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1 Capitolo 1 


Introduzione 


I beni architettonici storici e culturali rappresentano una risorsa inestimabile per il nostro Paese e la 
loro conservazione sta diventando sempre più una questione da affrontare, visto che sono soggetti 
allo scorrere del tempo e soprattutto ai continui eventi sismici che colpiscono le nostre regioni. 

Basti pensare al terremoto Umbria-Marche nel 1997, a L’ Aquila nel 2009, all’Emilia nel 2012 e al 
Centro Italia nel 2016, che hanno confermato soprattutto la grande problematica della vulnerabilità 
sismica delle chiese storiche in muratura. Il nostro Paese è un territorio sismico, con il quale si deve 
imparare a convivere, in modo da preservare il patrimonio culturale per le generazioni future, 
affrontando la necessità della riduzione del rischio sismico, promuovendo una cultura della 
conoscenza e della prevenzione. Al fine di sottolineare l’importanza del concetto di conservazione si 
può fare riferimento alle “Linee Guida per la valutazione e riduzione del rischio sismico del 
patrimonio culturale”, dove sono indicate tutte le procedure per intervenire con rispetto del bene 
culturale. 

Le chiese in muratura costituiscono un insieme estremamente vario e complesso per le tipologie e 
tecniche costruttive, per cui l’analisi del loro comportamento strutturale e la valutazione della loro 
sicurezza sono condizionate da notevoli incertezze nella definizione delle proprietà meccaniche dei 
materiali e delle condizioni di vincolo tra gli elementi. L’elevata vulnerabilità di questa tipologia di 
struttura nasce proprio dalla loro configurazione strutturale e geometrica, dai materiali eterogenei e 
dalla loro degradazione. In particolare, sono caratterizzate da alte pareti esterne molto larghe, senza 
pareti interne in direzione ortogonale, dove l’ampio spazio è spesso coperto da volte leggere e archi 
spingenti. 

Negli ultimi anni si sta affermando una tecnica in grado di valutare la risposta sismica e di controllare 
la funzionalità di una struttura nelle condizioni di esercizio, per poter progettare eventuali interventi 
di miglioramento sismico e rinforzo. Il monitoraggio strutturale dinamico (Vibration Based Structural 
Health Monitoring), a partire da una identificazione preliminare, permette di tenere sotto controllo 
l’evoluzione della dinamica della struttura. I dati forniti in tempo reale e in modo continuo sono un 
utile supporto ai fini della conservazione del manufatto e dell’ottimizzazione nella programmazione 
degli interventi di manutenzione. Inoltre, il monitoraggio consente di valutare tempestivamente la 
presenza del danno e può rivelarsi cruciale nel supporto del processo decisionale delle 
amministrazioni a seguito di eventi eccezionali. Infatti, il monitoraggio dinamico permette di 


controllare che il comportamento della struttura rimanga entro un determinato range “sano” e nel 


momento in cui ci sono delle anomalie, ovvero nel momento in cui vi è una variazione dei parametri 
dinamici, può significare che vi sono delle modifiche della struttura causate dalla presenza di un 
danneggiamento, poiché il comportamento dinamico dipende proprio dalle sue caratteristiche 
intrinseche (masse, rigidezze, gradi di vincolo ecc.). 

È proprio per questo che è uno strumento molto importante per le strutture storiche in muratura, 
poiché può costituire un campanello di allarme che permette di intervenire il prima possibile. Inoltre, 
la muratura storica è un materiale che presenta diverse incertezze, e nonostante si possano avere a 
disposizione i risultati di prove di caratterizzazione dei materiali, l’identificazione dinamica fornita 
dal monitoraggio consente di identificare la risposta globale della struttura, difficilmente ottenibile 
da indagini puntuali. 

Il limite del monitoraggio è che, avendo generalmente a disposizione un limitato numero di sensori, 
consente difficilmente la localizzazione del danno, ed è per questo che il più delle volte viene 
sviluppato un modello agli elementi finiti a supporto. Questo deve essere calibrato sulla base dei 
risultati di un’identificazione dinamica preliminare, e aggiornato in base ai risultati sperimentali 
provenienti dal monitoraggio, in modo che sia in grado di rappresentare al meglio il comportamento 
dinamico della struttura nel tempo. Una volta eseguito l’updating, il modello sarà in grado di fornire 
i dati mancanti per poter rilevare la presenza e la posizione del danno. In generale, il monitoraggio 
consente di definire nel tempo un pattern statistico rappresentativo delle condizioni operative della 
struttura; affiancando un modello a supporto del monitoraggio, la definizione del pattern statistico si 
definisce supervisionata. 

Il lavoro della presente tesi consiste nella creazione di un modello agli elementi finiti della chiesa di 
Santa Maria in Via a Camerino, gravemente danneggiata a seguito del sisma del 26 ottobre 2016, a 
supporto del sistema di monitoraggio in corso di installazione. Dopo avere eseguito un attento rilievo 
geometrico della struttura e le indagini visive necessarie per individuare le tipologie di murature 
presenti, è stato possibile disegnare un preciso modello tridimensionale, dal quale è stato ottenuto un 
modello agli elementi finiti solidi con il quale è stata poi eseguita l’analisi modale per poter 
individuare il suo comportamento dinamico. Grazie allo svolgimento di un’indagine dinamica 
preliminare, tramite prove di vibrazioni ambientali, è stato possibile avere dei dati sperimentali con i 
quali confrontare i risultati provenienti dal modello numerico, ed eseguire una calibrazione in modo 
da renderlo il più vicino possibile al comportamento reale. In particolare, l’attenzione è stata posta 
sull’ottimizzazione dei valori dei moduli elastici, sulla valutazione dell’effetto della rigidezza dato 


dalla copertura al tiburio e sul grado di collaborazione delle solette dei solai. 


2. Capitolo2 


Stato dell’arte 


2.1 Modellazione agli Elementi Finiti 


La modellazione agli elementi finiti, definita anche FEM (Finite Element Method), è una strategia di 
modellazione numerica adatta a qualsiasi tipo di struttura, largamente usata nell’ingegneria civile per 
capire e studiare il comportamento statico e dinamico delle strutture, anche con differenti tipi di 
materiali strutturali come il cemento, l’acciaio, il legno e la muratura. Essendo uno degli approcci più 
accurati trova un grande utilizzo nell’indagine di casi complessi, la cui valutazione in laboratorio 
porterebbe a costi eccessivi. 

Il suo obiettivo è quello di ottenere delle soluzioni approssimate e per far ciò il dominio della struttura 
viene suddiviso in tanti sottodomini, definiti elementi finiti, riducendo così il problema ad un numero 
finito di incognite. Discretizzando la struttura in un certo numero di elementi finiti, infatti, si definisce 
il campo incognito con una serie di funzioni approssimate, dette funzioni di forma, ognuna all’interno 
di ogni elemento. I singoli elementi possono avere forme diverse (triangolari, quadrilateri, esaedri e 
tetraedri) e possono essere monodimensionali, bidimensionali o tridimensionali. Sono collegati tra 
loro da dei nodi, detti punti nodali. Gli spostamenti di questi ultimi corrispondono ai gradi di libertà 
della struttura e su di essi si concentrano anche le forze esterne. 

All’aumentare del numero degli elementi finiti che compongono il dominio, dipendente dalla mesh 
di suddivisione utilizzata, il modello diventa più preciso e quindi il livello di approssimazione è 
condizionato dalle scelte di impostazione del modello. È importante dunque scegliere la giusta scala 
di modellazione, che nel caso specifico della muratura, può anche essere indipendente dalla reale 
tessitura. Le proprietà meccaniche dell’insieme giunti ed unità, infatti, possono essere rappresentate 
come un continuo omogeneo caratterizzato da un particolare legame costitutivo, consentendo una 
semplificazione del processo di discretizzazione. Inoltre, è possibile effettuare differenti livelli 
d’indagine in base a dove sia necessario approfondire in modo più accurato lo stato di tensione e 
deformazione, attraverso una suddivisione più fitta della mesh. 

La modellazione agli elementi finiti ha avuto negli ultimi anni sempre maggior diffusione uscendo 
dal campo della ricerca. In commercio si trovano diversi codici di calcolo che hanno reso questa 
strategia di modellazione sempre più reperibile e accessibile. In particolare, in questa tesi è stato 


utilizzato il software ANSYS. 


2.1.1 Scala di modellazione 


La muratura è composta da un insieme di blocchi congiunti tra loro con uno strato di malta o a secco. 
I blocchi (unità) e la malta sono due elementi differenti, ognuno con le sue proprietà meccaniche e 
data la discontinuità tra gli elementi, l’interazione tra i due è rappresentata da dei piani deboli definiti 
“giunti di interfaccia”. Questo comporta che il materiale è governato da una serie di meccanismi di 
rottura complessi che interessano i blocchi (Fig.1a) e la malta (Fig.1b, c), ma anche una combinazione 


tra i due (Fig.l d, e). 
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Figura 1: Meccanismi di rottura (Kumar et al., 2014) 


La precisione con cui si può riuscire a cogliere il reale meccanismo di rottura della muratura è quindi 
dipendente dalla tipologia di discretizzazione, che può essere suddivisa in tre categorie: 

- Micro-modellazione; 

-  Micro-modellazione semplificata; 

- Macro-modellazione. 
Nella micro-modellazione vengono rappresentati i singoli componenti della muratura, ognuno con le 
proprie proprietà meccaniche, rappresentando in modo molto dettagliato la muratura, ovvero la sua 
reale tessitura. Con questa strategia si riescono ad individuare con una certa precisione tutti i 
meccanismi di rottura che interessano i blocchi e la malta, ponendo particolare attenzione alla 
modellazione dei giunti, poiché si osservato che la superficie di contatto è irregolare e concentrata 
nella parte centrale, probabilmente a causa del ritiro della malta e della tecnica di posa. 
Essendo una tipologia molto dettagliata, e considerando la grande complessità degli edifici in 


muratura, non risulta sempre applicabile questo tipo di modellazione. Tuttavia, potrebbe essere 


un’ottima soluzione per piccole porzioni murarie, sempre con tessitura abbastanza regolare, dove è 
necessario porre l’attenzione e fare specifiche analisi. 

La micro-modellazione (Fig.2) si divide in dettagliata, appena descritta, e semplificata, chiamata 
anche modellazione ad interfacce, nella quale vengono trascurati i giunti di malta e si prendono in 
considerazione solo i blocchi della muratura collegati tra loro dalle superfici di interfaccia. Dato che 
in questo modo lo spessore dei giunti di malta viene trascurato, le dimensioni delle unità vengono 
scalate. È quindi una semplificazione che permette di essere applicata a più tipologie di muratura, 


considerando però una leggera perdita di accuratezza dei risultati. 


Mortar Unit “ Unit” 
Interface 
UnivMortar 


“Joint” 


Figura 2:Micro-modellazione a sinistra, micro-modellazione semplificata a destra, 
La macro-modellazione (Fig.3), chiamata anche omogenea, considera invece la muratura come un 


unico materiale anisotropo, trascurando completamente la distinzione tra unità, giunti di malta e 


Fa 


superfici di interfaccia. Composite 


Figura 3: Macro-modellazione, Laurenco, 2002. 


In questo modo, si raggiunge un minore livello di accuratezza, poiché la mesh è indipendente dalla 
reale tessitura muraria. Tuttavia, è una metodologia molto utilizzata nell’ambito ingegneristico, 
poiché quando ci si riferisce ad un’intera struttura, è necessario cogliere il danneggiamento su grande 
scala, trascurando la precisa localizzazione degli sforzi. Inoltre, la maggior parte di strutture storiche 
in muratura presenta una tessitura irregolare, caratterizzata dalla presenza di blocchi disomogenei e 
molto spesso, dalla presenza di murature a sacco. Questa tipologia di muratura è di difficile 
valutazione, e quindi è improbabile pensare di modellarla con una micro-modellazione. 

Il caso studio analizzato in questa tesi è molto complesso, sia per la geometria che per la tipologia di 


muratura, ed infatti è stato modellato tramite una macro-modellazione. 


Quindi, una strategia di modellazione non può essere preferita ad un’altra, poiché esistono differenti 
campi di applicazione per la macro e la micro-modellazione. 

Ci sono degli studi 1, basati sulla comparazione di esperimenti eseguiti in laboratorio, che dimostrano 
l'efficacia della macro-modellazione nel rappresentare il comportamento della muratura, a volte 
maggiore della micro-modellazione, poiché anche se più dettagliata, presuppone la gestione di un 
maggior numero di parametri meccanici che possono influire notevolmente il risultato. Il modello 
omogeneo invece necessita la determinazione di proprietà medie che rispecchiano il comportamento 
macroscopico del solido murario, anche se di non facile determinazione. 

La stima delle caratteristiche meccaniche elastiche macroscopiche della muratura si basa sulla 
conoscenza delle proprietà elastiche della malta e dei blocchi, attraverso tecniche di 
omogeneizzazione 0 osservazioni fenomenologiche. In questo studio, dato il basso livello di 
conoscenza dei materiali, si è partiti dall’individuazione della tipologia di muratura presente e si è 


considerato il valore medio delle caratteristiche elastiche delle tipologie presenti in letteratura. 


2.2. Model Updating con prove di vibrazione ambientale 


Il Model Updating di un modello agli elementi finiti è una procedura molto efficace per determinare 
i parametri sconosciuti o incerti di un modello strutturale ed ha lo scopo di minimizzare le differenze 
tra il modello sperimentale e il modello numerico. 

Le strutture in muratura sono caratterizzate, la maggior parte delle volte, dall’avere un’importante 
incertezza nei confronti dei materiali, delle tipologie costruttive, delle masse presenti e delle 
condizioni al contorno, che rendono difficile la costruzione di un adeguato modello numerico che 
rappresenti fedelmente il loro comportamento statico e dinamico. È proprio per questo che i risultati 
numerici dovrebbero essere validati e confrontati con le misure sperimentali. Negli ultimi anni questa 
procedura di updating sta ricevendo sempre più attenzione e significato poiché i modelli aggiornati 
possono essere usati per la localizzazione del danno e a supporto di un monitoraggio della salute 
strutturale. 

Le analisi modali sperimentali forniscono le caratteristiche dinamiche delle strutture, come le 
frequenze naturali, i modi di vibrare e i fattori di smorzamento, che forniscono un meccanismo di 
controllo per validare i risultati provenienti dal modello sperimentale agli elementi finiti, poiché sono 
prove eseguite su strutture reali in sito e quindi 1 risultati possono essere accettati come 1 principali 
parametri di riferimento. 


Sono due le principali tecniche con cui si possono eseguire le analisi modali sperimentali: 


- FVT (Forced Vibration Test), ovvero prove forzate di vibrazione basate su una forza di input 
nota, usando un’eccitazione artificiale come un impulso hammer o uno shaker hydraulic; 
-  AVT (Ambient Vibration Test), ovvero prove con vibrazioni ambientali basate su una forza 


di input sconosciuta, come il vento, il traffico e il terremoto. 


L’aggiornamento del modello agli elementi finiti può essere eseguito secondo due metodi: updating 
del modello manuale e l’updating del modello automatico. 

L'aggiornamento del modello manuale consiste nel cambiare manualmente la geometria e i parametri 
eseguendo diverse prove e valutando di volta in volta l’operazione da eseguire con giudizio 
ingegneristico. La parte più importante di questa procedura consiste nell’individuazione dei parametri 


incerti di aggiornamento, ovvero: 


- Proprietà dei materiali: modulo di Young, peso specifico, coefficiente di smorzamento; 
- Proprietà geometriche: sezioni e masse aggiuntive; 
- Proprietà di modellazione: condizioni al contorno. Solitamente le più usate come parametri di 


aggiornamento sono le proprietà dei materiali e le condizioni al contorno. In particolare, è il 
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modulo di Young il parametro più utilizzato proprio perché è il più incerto delle strutture 


storiche in muratura. 


La procedura di aggiornamento automatico considera molti più parametri incerti in modo da 
ottimizzare ancora di più la correlazione tra 1 due modelli e consiste solitamente in strategie di 
aggiornamento globali e locali. La prima considera un singolo valore per ogni parametro in tutti i 
modelli, e viene fermata quando viene raggiunta una determinata tolleranza. La seconda invece, 
considera le proprietà del materiale a differenti valori per ogni elemento di mesh e viene eseguita 
dopo l’aggiornamento del modello globale per eseguire analisi specifiche. 
I principali step, schematizzati anche in figura 4, sono: 

- Selezione dei parametri incerti; 

- Determinazione dell’efficacia di ogni parametro; 

- Specifica dei limiti ammissibili per ogni parametro; 


- Studi di convergenza per minimizzare 1 risultati sperimentali e numerici. 


Finite Element Model Ambient Vibration Testing 


Manual Model Updating 


Global Updating 


Selection of Uncertain 
Parameters 


Changmng of Parameters 


by Trial and Error 


Bayesian Parameter Estimation 


YES 


Figura 4: Esempio di un Flowchart del model updating del bastione Zagonos in Turchia. 
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2.2.1 Modal Assurance Criterion (MAC) 


Per valutare e avere un’indicazione sulla distanza dei vari modi di vibrare ottenuti sperimentalmente 
e numericamente, viene usato il Modal Assurance Criterion (MAC). È un buon indicatore statistico, 
poiché è più sensibile alle grandi differenze e insensibile alle piccole differenze tra 1 modi di vibrare, 
ovvero riesce a valutare anche modi di vibrare molto vicini e a dare una buona valutazione sulla loro 
corrispondenza. 

Viene calcolato come un prodotto scalare tra due insiemi di vettori {Q@A} e {@58}. Gli scalari risultanti 
sono disposti nella MAC matrix. 

2 2 
acta) = — PATRA 
{0,3} (0 Da) 

Assume valori compresi tra 0 e 1, dove 1 indica la piena corrispondenza dei modi mentre un valore 
vicino lo 0 indica che non vi è corrispondenza. 
Può assumere valori prossimi allo zero per i seguenti motivi: 

- Il sistema non è stazionario a causa di cambiamenti nella massa, rigidità e proprietà di 

smorzamento durante il test; 

- Il sistema non è lineare; 

-  C’è un disturbo sul modo di vibrare di riferimento; 

- La tecnica di estrazione dei parametri non è valida per l’insieme di dati misurati; 

- Le forme modali sono linearmente indipendenti. 


Invece, assume valori prossimi all’unità quando: 


- Il numero dei gradi di libertà di risposta non è sufficiente per distinguere le forme modali 
indipendenti; 

- Le forme modali sono il risultato di forze non misurate sul sistema; 

- Le forme modali sono principalmente disturbo coerente; 


- Le forme modali rappresentano lo stesso movimento, differente solo da uno scalare. 


Il MAC può indicare solo una corrispondenza, non la validità di un modello, così è principalmente 
usato per una valutazione preliminare. Infatti, in questa tesi viene usato questo indicatore per avere 
un’idea della corrispondenza tra i modi sperimentali e i modi numerici ad ogni step di aggiornamento 


del modello. In figura 5, sono mostrate le varie tipologie di rappresentazione della MAC. 
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Figura 5: Rappresentazione numerica, 2D e 3D dei valori MAC. 
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0.025 
0.023 
0.067 


2.3. Casi studio 


I casi studio analizzati si dividono in cinque casi relativi alla modellazione agli elementi finiti e cinque 
casi inerenti l’updating di un modello FEM. La localizzazione dei casi studio è mostrata nelle figure 


6 e7. 


11-59 quilà 
4 1 3 i a 
UE (1.48 113 


2: 


Figura 6: Localizzazione casi studio, Europa 


2 + 


Figura 7: Localizzazione casi studio, America del Sud 
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2.3.1 Casi studio FEM 


2.3.1.1 Abbazia di Farneta 


STATO Italia 
CITTÀ Val di Chiana (AR) 
ANNO DI COSTRUZIONE 700-800 d.C. 


https://www.researchgate.net/p 
ublication/223444850 


FONTE 


Tabella 1: Identi i hiesa. 
cita RARE RA Figura 8: Vista dall'alto, Google Maps. 


Figura 9: Vista frontale. Figura 10: Vista posteriore. 


Cenni storici e descrizione della struttura 


La struttura originale fu costruita tra il 700 e 1° 800 d.C. dal conte di Ronzano. Raggiunse il periodo 
di massima maestà tra il decimo e il tredicesimo secolo, quando fu occupata dai monaci benedettini. 
Dopo questo tempo, l’ Abbazia perse la propria autonomia, sia civica che ecclesiastica, fino agli inizi 
del diciottesimo secolo. Oggi è sotto la giurisdizione ecclesiastica di Cortona. 

Nel corso degli anni la struttura originale romanica ha subito diverse aggiunte e modifiche, infatti 
l'Abbazia attuale presenta delle differenze rispetto all’originale, come ad esempio la lunghezza ridotta 
della navata principale. Questa riduzione è stata fatta durante il diciottesimo secolo quando, 
probabilmente, la costruzione subì un crollo e quindi fu necessaria una ristrutturazione. Un altro 
esempio che testimonia le varie modifiche fatte nel corso degli anni, è la presenza di diverse tipologie 


99 66 


di tessitura muraria, differenti per materiali (mattoni, pietre, etc.) e per forma (‘opus incertum”, “opus 
mixtum”, “opus quadratum”, etc.). 
Nel 1964 la chiesa subì una restaurazione in cui venne costruita una nuova fondazione, andando così 


a modificare il tenore idraulico del terreno, e in cui venne svuotata la cripta, piena di materiale 
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accumulato negli anni. Questi lavori portarono alla 
formazione di diverse lesioni nella parete settentrionale 
del transetto, in prossimità dell’abside. 

L’abbazia di Farneta ha la planimetria tipica di 
un’ampia classe di chiese romane esistenti nel centro 
Italia (figura 11 e 12), poiché presenta una navata 
principale, un cleristorio e una tripla abside. Le 
dimensioni principali consistono in una lunghezza 
massima di circa 26,5 m, una larghezza massima di 


circa 21,5 m e un’altezza di circa 11 m. Lo spessore 


delle pareti in muratura varia tra 0,70 m (pareti della 


Figura 12:Planimetria della chiesa. 


HYS 


navata principale) e 0,90 m (cleristorio). Dal lato 
sinistro e destro del cleristorio è possibile scendere 


tramite delle scale, in una cripta interamente sviluppata 


sotto la zona del cleristorio. Il tetto della chiesa è in 
legno, mentre nella cripta sono presenti delle volte, 
composte da mattoni e pietre, connesse con i muri di Figura 11:Planimetria della cripta. 

confine e sostenute da due colonne circolari in granito. Nonostante la presenza di differenti tessiture 
murarie, la costruzione è prevalentemente in muratura di arenaria irregolare, materiale che si trova in 
zona, con giunti di malta di calce spessa. Le proprietà meccaniche dei materiali sono state ricavate da 


dati sperimentali di strutture simili e da prove in situ. 


Modellazione FEM e analisi svolte 

La chiesa in esame viene analizzata per valutare il comportamento strutturale e la vulnerabilità 
sismica allo stato di conservazione attuale, e inoltre viene analizzato il comportamento rispetto a 
diverse proposte di rafforzamento, tramite il metodo agli elementi finiti. Vengono realizzati due 
modelli per le analisi. Il primo descrive la struttura allo stato attuale, “as is” (Fig. 13), in cui non 
vengono considerati i collegamenti tra le varie strutture murarie e il tetto poiché non sono efficienti 
nella realtà. Nel secondo vengono considerati i collegamenti, visto che rappresenta la struttura dopo 


ipotetici interventi di rafforzamento. 
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Software Ansys 


Tipologia di Elemento Elemento solido tetraedrico con 8 nodi 
*Solid45 
eSolid65 


Dimensione *10620 nodi 

*5963 elementi 

Utilizzo *Valutazione comportamento strutturale e vulnerabilità 
sismica; 

*Valutazione efficacia interventi; 


Tabella 2: Descrizione modello FEM. 


1 ANSYS 5.5.2 


ELEMENT SOLUTION 
STEP=1 

SUB =1 

TINE-1 

sY (NORVG) 
RSYS=0 
PoverGraphies 
EFACET=1 


Farneta 


Figura 13:Sforzo verticale, modello "as is". 


Si è partiti dall’analisi statica lineare, in cui le pareti sono state modellate mediante elementi solid45 
(elementi lineare elastici isoparametrici a 8 nodi) considerando il materiale omogeneo e isotropo con 
un comportamento lineare elastico. In questa fase la struttura è soggetta solamente ai carichi verticali, 
ovvero il peso proprio e i carichi applicati in sommità delle pareti provenienti dal tetto. 
Successivamente si è passati all’analisi statica non lineare per indagare sull’origine dell’attuale stato 
fessurativo. In questo caso gli elementi utilizzati nel modello sono i solid65 e il comportamento del 
materiale viene descritto dal criterio di rottura Drucker-Prager perfettamente plastico. 

Per valutare il comportamento sotto l’azione sismica è stata effettuata un’analisi sismica, dove il 


modello è stato modificato aggiungendo degli elementi contact52, intorno all’area dove sono presenti 


Figura 14:Primi tre modi di vibrare del modello "as is". 
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le lesioni e per la combinazione dei carichi è stata usata la tecnica SRSS (Square Root of Sum 


Squares). Nelle figure 14 e 15 vengono riportati i risultati dell’analisi dinamica. 


DISPLACEMENT pe” mf 


Figura 15: Primi tre modi di vibrare del modello usato nell’analisi sismica 


Ansrs 5,5.2 
Farneta abbey (“as is”) Farneta abbey (nonlinear analysis) cRACEI AND cRUSEID 

Mode n° —Period(s) Period (5) Bea Tn. 
e, cri... rrrrr1r’’’’’. poca Ta 

1 0.318 0.479 î | 

2 0.305 0.428 | 1 

3 0.273 0.350 \ 

4 0.225 0.315 

È) 0.181 0.292 | 

6 0.169 0.255 L 

7 0.169 0.252 

8 0.162 0.246 

9 0.160 0.218 

10 0.141 0.203 

Tabella 3: Confronto dei principali periodi del modello "as Figura 16:Percorso fessurativo al collasso per 

is" e modello usato nell'analisi sismica. il modello rinforzato (A). 


Infine, sono state svolte nuovamente delle analisi non lineari, però incentrate sul secondo modello, 
ovvero quello rinforzato. Un esempio viene mostrato in figura 16. 


Sono state progettate tre ipotesi di intervento: 


e Modello (A): connessione tra copertura e pareti; inserimento di catene; nuova fondazione. 
e Modello (B): rafforzamento modello (A). 
e Modello (C): rafforzamento modello (A); inserimento cordolo in fondazione lungo 1 lati 


dell’edificio. 


b) 


|rarneta rarneta 2 
| I E OS 


Figura 17: a) Modello rinforzato (A), b) Modello rinforzato (C). 
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Risultati e conclusioni 


Le analisi del comportamento statico e sismico 


lm Identified building MI Strengthening [A] Strengthening [B] O Strengthening [IC] | 


della chiesa hanno mostrato due principali 


debolezze della struttura. La prima è la 


mancanza di connessioni tra i muri ortogonali, 


Period (s) 


mentre la seconda è la necessità di consolidare 


il terreno intorno l’abside. Così si è scelto di 


proseguire le analisi per andare a valutare i 
Mode Shape 


benefici dati dagli interventi di rafforzamento hail 
(figura 18). Queste analisi sono state svolte grazie alla costruzione di un modello agli elementi finiti, 
grazie al quale si possono ottenere informazioni significative rispetto alla comprensione dei danni 


esistenti e alla progettazione minima e adeguata di rafforzamento. 


2.3.1.2 Complesso di San Francesco 


STATO Italia 
CITTÀ Cagli (PU) 
ANNO DI COSTRUZIONE 1234 d.C. 


https://www.researchgate.net/pub 
lication/308153137 


FONTE 


E is MIT ; 
7 SL! 7 ®- 


Figura 20: Vista dall'alto, Google Maps. 


Tabella 4: Identificazione chiesa. 


Figura 19: Complesso di San Francesco 


Cenni storici e descrizione della struttura 


La conoscenza dell’evoluzione storica dell’edificio è fondamentale per la determinazione delle 


interazioni tra le parti che compongono l’intero complesso. 
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Il convento, la chiesa e il campanile vennero costruiti nel 


lo] 1234 - existing original core 


1891 - first addition 


1234 e l’intero complesso venne devastato da un incendio 


demolished old sacristy 


XX century - new constructions 
(structurally independent) 


E 
nel 1415 e ricostruito 5 anni dopo, con l’aggiunta di una È 
seconda elevazione in alcune aree dell’edificio (Fig. 21). È 

La chiesa è caratterizzata da uno sviluppo longitudinale e 
interagisce con la struttura del monastero, il campanile e 


l’area dell’abside. L'intervento strutturale più interessante 


VIA SAN FRANCESCO 


da ricordare è stato approvato nel 1891 che portò 
all’aggiunta di due bracci per formare una struttura 
simmetrica, ovvero l’attuale configurazione. 

n Figura 21: Configurazione del monastero nei secoli 
E stato eseguito un dettagliato rilievo geometrico e attente 

ispezioni visive per ricavare informazioni necessarie alla costruzione di modelli numerici (Numerical 

Models-NMs) usati per valutare il comportamento statico e sismico della struttura. Inoltre, è presente 

una valutazione delle caratteristiche meccaniche dei materiali (Tab.5), la cui incertezza è stata presa 


in considerazione tramite l’utilizzo del fattore di confidenza (FC) e del livello di conoscenza (LC), 


come dettato dalla normativa (NTC18 e Eurocodice 8). 


Figura 22:Due prospetti del monastero. 


fr fi Y E v Ge Gr Confidence factor (CF) 
[MPa] [MPa] [N/mm}] [MPa] [N/mm] [N/mm] 
Stones split with good weaving masonry - Material 1 
KLI 3.255 0.3255 0.000021 1740 0.3 2.6436 0.026436 1.35 
KL2 4.507 0.4507 0.000021 1740 0.3 3.3199 0.033199 120 
KL3 6.422 0.6422 0.000021 1740 0,3 4.2540 0.04254 1.00 
Squared block stone masonry - Material 2 
KL1 5.333 0.5333 0.000022 2800 0.3 3.7353 0.037353 1.35 
KI 7.000 0.7000 0.000022 2800 0.3 4.5185 0.045185 1.20 
KL3 9.600 0.9600 0.000022 2800 0.3 5.6366 0.056366 1.00 
Solid bricks and lime mortar masonry - Material 3 

KLI 2.667 0.2667 0.000018 1500 0.3 22994 0.022994 1.35 
KL2 4.000 0.4000 0.000018 1500 0.3 3.0540 0.030540 120 
KL3 6.000 0.6000 0.000018 1500 0.3 4.0563 0.040563 100 


Tabella 5: Caratteristiche meccaniche dei principali elementi 
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Modellazione FEM e analisi svolte 

Il comportamento globale della struttura e le proprietà dinamiche del complesso sono stati valutati 
tramite una modellazione agli elementi finiti, dove il comportamento non lineare della muratura è 
stato preso in considerazione dall’assunzione di specifici legami costitutivi. In particolare, l’obiettivo 


principale dello studio è quello di valutare l’influenza, nella risposta globale della struttura, da: 


e lascelta del punto di controllo; 
e larigidezza dei piani (effetto diaframma); 


e illivello di conoscenza, così definito dalla norma. 


L’ultima analisi consiste nel confronto tra i risultati ottenuti con il 3D FEM e il semplificato 3D EFM 


(Equivalent Frame Model). 


Software Midas FEA 
Tipologia di Elemento Elemento solido tetraedrico con 4 nodi 
Dimensione *Modello 1 — 79832 nodi, 268328 elementi, 231549 dof. 


«Modello 2 > 79841 nodi, 268328 elementi, 229773 dof; 
«Modello 3 —> 79914 nodi, 268328 elementi, 221772 dof. 


Utilizzo *Valutazione comportamento strutturale e vulnerabilità 
sismica; 

*Confronto con altra tipologia di modellazione (EFM). 
Tabella 6: Descrizione modello FEM. 


Le pareti in muratura sono modellate con un elemento solido tetraedrico con 4 nodi e per la 
discretizzazione è usata una mesh regolare ottimizzata. Il comportamento non lineare è rappresentato 
da Total Strain Crack Model (Fig. 23) disponibile nel codice di calcolo utilizzato, e in figura possiamo 
vedere i legami costitutivi usati. Il modello riproduce accuratamente la geometria delle strutture, 
focalizzandosi sulle variazioni dello spessore dei muri, sulle irregolarità delle strutture geometriche e 


strutturali e sulle connessioni dei muri. 


a) b) 


Figura 23:Leggi costitutive tensione-deformazione: a) 
compressione uniassiale della muratura; b) tensione uniassiale 
della muratura. 
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Figura 24:Modello agli elementi finiti di San Francesco in Cagli. 


Per analizzare l’effetto della rigidezza dei 


diaframmi sono stati realizzati tre modelli diversi: 


e Modello 1: solai deformabili (a); 
e Modello 2: solai rigidi solo in una parte, 
dove è presente la scuola (b); 


e Modello 3: solai rigidi (c). 


Viene svolta un’analisi modale e nella figura 25 vengono riportati i principali modi di vibrare dei tre 
modelli. Si nota che la prima frequenza non è associata allo stesso elemento del complesso e significa 
che la parte più vulnerabile del monastero cambia con l’inserimento dei diaframmi rigidi e con il 


miglioramento del comportamento scatolare della struttura. 
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Mode 2 Mode 18 Mode 26 Mode 29 Mode 40 
T=0.788 T=0.397 T=0.256 T=0.224 T=0.202 


Mode 25 cia 33 Mode 34 Mode 38 Mode 49 
T=0.255 T=0.199 T=0.189 =0. T=0.171 T=0.149 
Mode 3 Mode 7 Mode 8 Mode 10 Mode 11 Mode 12 
T=0.300 T=0.192 T=0.183 T=0.163 T=0.153 =0.147 


Figura 25:Risultati delle analisi modali per i tre modelli: a) Model 1, b) Model 2, c) Model 3. 


Vengono svolte anche analisi statiche non lineari, dalle quali si ricavano le curve di capacità in base 
ai carichi nelle diverse direzioni. Vengono predisposti gli stessi punti di controlli in tutti 1 modelli, 
così è possibile valutare la variazione della risposta sismica cambiando la rigidità dei solai e la 


capacità di resistenza. 


Per eseguire analisi globali con un onere computazionale ragionevole e pochi parametri meccanici, 
viene realizzato anche un modello EFM, che consiste nel rappresentare le strutture come un 
assemblaggio di elementi strutturali macroscopici. In questo caso il software utilizzato è 3MURI. I 


due modelli sono stati messi a confronto. 


Risultati e conclusioni 


I risultati di questo lavoro confermano le vulnerabilità tipiche di questo tipo di edifici, come la 
posizione eccentrica del campanile, piani flessibili, l’assenza di pareti trasversali nella chiesa e la 
bassa capacità dissipativa delle colonne del cortile, che possono influenzare il comportamento sismico 
globale di questo edificio. 

Grazie a una buona conoscenza del manufatto è stato possibile costruire un accurato modello agli 
elementi finiti, non molto semplice da ottenere vista la grande complessità della struttura. Questo ha 
garantito un ottimo livello di accuratezza dei risultati insieme a un buon grado di dettaglio relativo 
allo stato fessurativo, anche se richiede tempi di calcolo molto lunghi. 

Infine, si è visto che i risultati ottenuti con 3MURI hanno una buona correlazione con i risultati 
estrapolati da Midas FEA e con tempi di calcolo minori. Tuttavia, per ottenere risultati più precisi 


bisogna utilizzare il FEM. 
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2.3.1.3 Chiesa di Santa Maria della Carità 


i r——_ 


fe GastronomicatS:!Martino li 
Di Dimarco Miriam 


STATO Italia 
CITTÀ Ascoli Piceno (AP) 
SODI Sconosciuto 
COSTRUZIONE 


https://www.researchgate.net/publicatio 
n/320208593 


FONTE 


Tabella 7: Identificazione chiesa. 


SIENA ss po 


Cenni storici e descrizione della struttura 


Il periodo di costruzione della chiesa non è noto, ma in una pergamena viene nominata una piccola 
chiesa costruita nel 1306. Nel corso del sedicesimo secolo l’edificio è stato parzialmente demolito e 
ricostruito, in cui si modificò la facciata con l’innalzamento del timpano e con l’apertura dei due 
portoni frontali. A seguito di questo intervento la navata centrale è stata coperta da una serie di volte 


in muratura, collegate da un arco trionfale con l’abside, costruito nello stesso periodo. Il campanile 


attuale fu costruito verso la fine del diciassettesimo secolo, in una posizione diversa dal precedente, 
| | 


come attestano dei documenti iconografici. | 


Figura 28: Pianta (sinistra) e sezione (destra) della chiesa 
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La chiesa ha un’unica navata di 12,93 m di larghezza, 23,15 m di lunghezza e 16,5 m di altezza. 
L’abside si trova dentro la sacrestia, che è divisa in due piani, con un’altezza massima di 15,90 m. 
Dalla sacrestia si accede al campanile, che ha una sezione quadrata di lato 4,31 m e un’altezza 
massima di 27,80 m. Le pareti hanno uno spessore che varia tra 2,21 m (pilastri della navata) e 0,91 
m (muri della navata). 

La chiesa presenta differenti tipologie di murature (mattoni, travertino, pietra etc.) e in particolare le 
pareti della navata sono caratterizzate dalla presenza di due differenti murature in elevazione, dovute 
a un differente periodo di costruzione: fino a 11,08 m è presente il travertino, mentre da lì in su pietra 
irregolare. 

Le caratteristiche meccaniche dei materiali sono state ricavate da prove in situ e da norme per edifici 


esistenti in muratura. 


Modellazione FEM e analisi svolte 


Lo scopo è quello di andare a valutare la vulnerabilità sismica della chiesa dopo la chiusura a causa 
del terremoto Marche-Umbria-Lazio-Abruzzo avvenuto nel 2016, tramite l’utilizzo di un modello 
agli elementi finiti (figura 30). Il modello FEM è stato costruito per riprodurre la geometria della 
struttura, concentrandosi sulle variazioni dello spessore delle pareti, sulle irregolarità geometriche e 
strutturali e sulle connessioni tra i muri. Si considera una completa separazione tra la chiesa e 


l’annesso, poiché è stato costruito in un’epoca successiva. 


MO] LO ] 
Software Midas FEA 


Tipologia di Elemento |Elemento solido tetraedrico con 4 nodi 
Dimensione *17136 nodi; 

e59632 elementi; 

*498227 dof. 
Utilizzo Valutazione vulnerabilità sismica. 


Tabella 8: Descrizione modello FEM 


I solai vengono considerati flessibili e il modello 
considera il contributo di un intervento eseguito nel 
2010, dopo il terremoto dell’ Aquila del 2009. Inoltre, per 


le analisi statiche non lineari, il modello è considerato 


incastrato alla base. 


a) b) 


Il comportamento non lineare delle pareti in muratura è Figura 29:Leggi costitutive tensione-deformazione: a) 
compressione uniassiale della muratura; b) tensione 


rappresentato dal Total strain crack model presente lo Re ariana. 
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Base shear (kN) 


all’interno del programma utilizzato. Il comportamento di trazione/compressione è stato modellato 


secondo 1 legami costituivi illustrati in figura 29. 


Figura 30: Modello agli elementi finiti della chiesa. 


Il comportamento sismico è stato analizzato con un’analisi pushover e le principali curve di capacità 
ottenute sono riportate in figura 31. Si può notare che nella direzione X la struttura ha un 
comportamento fragile, mentre nella direzione Y si osserva un comportamento più duttile. Nella 
figura 32 vengono riportati 1 cracking pattern nelle due diverse direzioni e si può vedere che in X 
l’evoluzione del danno è concentrata intorno alla torre, vicino alla connessione tra l’arco trionfale e 
le pareti della navata, mentre in Y si notano danni più diffusi nella volta principale e nella sua 


connessione con le pareti della navata. 


55% Fb"... 834.97KN 47% Fb max 390, 26KN 
Fagade control point 
——= PushMass_-X 
= — = PushMass_-Y 
100% Fd", 1512.84kN 100% Fb*,na,830.35kN 
0.004 0.006 0.008 0.01 
Displacement (m) 
Figura 31: Curve di capacità. 


-20% Fb" 1263.75kN 


-6% Fb* nm 783. 19kN 


a) 
Figura 32: Cracking pattern ottenuti dalle analisi Pushover: 


carico uniforme in direzione X (a), e in direzione Y (b). 


b) 
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Risultati e conclusioni 


Le analisi sismiche consentono di avere un quadro degli elementi più vulnerabili della struttura, grazie 
anche all’utilizzo di una modellazione agli elementi finiti. In questo caso sono l’eccentricità della 
torre e la volta principale. Queste informazioni possono essere utilizzate per progettare lavori di 
adeguamento locale e globale, oltre ad essere di riferimento per lo studio dell’ampia varietà di chiese 


storiche presenti. 


2.3.1.4 Chiesa di San Pietro di Coppito 


STATO Italia 
CITTÀ L'Aquila (AQ) 
ANNO DI COSTRUZIONE XII-XIII secolo 
https://usra.it/wp- 
FONTE content/uploads/2019/10/252- 
1347-2206-1-DR.pdf 


Tabella 9: Identificazione chiesa. 


Figura 34:Chiesa di San Pietro di Coppito 


Cenni storici e descrizione della struttura 


La costruzione della chiesa risale alla metà del secolo 
XIII e nel corso degli anni subisce diverse 
trasformazioni, la prima delle quali dopo il terremoto 
del 1315: di questo periodo sono ad esempio i tre 
absidi, il prospetto di sinistra e probabilmente il 


portale principale. In seguito, fu sottoposta ad altri due 


ampliamenti, tra cui la realizzazione di una cappella 05 15 30m 
Figura 35: Planimetria della chiesa. 
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alla destra del transetto e forse una seconda cappella, di cui non è rimasta documentazione. Dopo il 
terremoto del 1703 subì gravi danni strutturali: il crollo della facciata e della torre campanaria e la 
distruzione delle due cappelle appena citate. La ricostruzione portò a una radicale modifica 
dell’impianto interno, con l’aula unica che venne allungata fino all’abside maggiore. Tra gli ultimi 
interventi alla chiesa ricordiamo quello avvenuto tra il 1969 e il 1971. 

Con il terremoto del 2009 fu nuovamente danneggiata, subendo altri crolli e danni strutturali. 

Dal rilievo geometrico si osserva una pianta asimmetrica (figura 35), con un’aula principale di 
dimensioni importanti, alla quale si affianca una piccola navata laterale, che termina poi in un 
transetto molto ampio dal quale si accede alle tre absidi a pianta poligonale. Da quella di destra si può 
accedere alla torre campanaria a base ottagonale. Il transetto è diviso in due parti da un arco trionfale 
costituito da tre arcate che s’impostano su due pilastri di sezione ottagonale. In elevazione la struttura 
risulta irregolare, generando una diversa distribuzione delle masse e delle rigidezze. La facciata non 
è vincolata negli angoli superiori, favorendo il ribaltamento fuori dal piano e il crollo parziale della 
parte sommitale della stessa, e riducendo l’efficacia effettiva del comportamento scatolare 
complessivo. 

La caratterizzazione meccanica dei materiali è stata resa possibile tramite prove effettuate su carote, 
prelevate in diverse parti della struttura, fornendo in questo modo dati fondamentali per la corretta 


costruzione di un modello agli elementi finiti. 
Modellazione FEM e analisi svolte 


L’oggetto dello studio è la valutazione della risposta sismica della chiesa a seguito dei danni riportati 
con il terremoto del 2009. Lo scopo dell’analisi è quindi quello di identificare i principali modi di 
vibrare a cui è soggetta la struttura tramite un modello agli elementi finiti, per poi individuare i 
principali meccanismi di collasso che interessano i diversi macro-elementi, che identificano a loro 
volta il complesso strutturale, al fine di confrontarsi poi con il reale danno subito. 

Vengono realizzati due modelli FEM, uno per l’intera struttura (figura 36) e l’altro per l’arco trionfale 


(figura 38). 


Software Modello intera struttura Modello arco trionfale 
Strauss 2.3 DIANA 9.4.2. 

Tipologia di Elemento |Brick Non specificato 

Utilizzo eAnalisi modale; Analisi del meccanismo di 


eConfronto con meccanismi | danno. 
di collasso dei 


macroelementi. 
Tabella 10: Descrizione modello FEM. 
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Figura 36: Modello agli elementi finiti chiesa di San Pietro di Coppito 


Si riportano, in forma di scheda, i confronti fra danno rilevato, meccanismo di crisi proposto dalle 


linee guida e i modi di vibrare derivanti dall’analisi agli elementi finiti. 


; 4 Meccanismo n°1 dell’abaco 
Lesione tra la facciata e l'aula - pica 
dei meccanismi di collasso 


Meccanismo di collasso dell’arco trionfale 


sno mm EER Meccanismo n°13 dell'abaco 
Lesione sull arco trionfale È PA 7 
dei meccanismi di collasso 


Modo 1 in direzione X. massa partecipante: 5,6% Modo 4 in direzione Y, massa partecipante: 10,2% 


o) Meccanismo della sommità della facciata Ribaltamento della torre campanaria 


(3 i 


Meccanismo n°27 dell’abaco 
Meccanismo n°2 dell'abaco dei meccanismi di collasso 
dei meccanismi di collasso ni 


Crollo della sommità della facciata 


Modo 22 in direzione X, massa partecipante: 6,78% ù Modo 17 in direzione Y. massa partecipante: 5.12% 
Figura 37: Schede di confronto: a) meccanismo di danno della facciata; b) Meccanismo di collasso dell'arco 
trionfale; c) Meccanismo di collasso della sommità della facciata; d) Meccanismo di collasso della torre campanaria 
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In seguito al sisma dell’aprile 2009, l’arco trionfale ha prevalentemente lavorato nel proprio piano 
subendo lesioni significative, ma senza collassare. Si è scelto di utilizzare un modello FEM, poiché 
un’analisi manuale del cinematismo risulta più complessa, anche per la presenza dell’arco rampante. 
Nel modello si tiene conto delle condizioni di vincolo che la struttura fornisce all’arco stesso, quindi 
viene impedita la traslazione fuori dal piano, in corrispondenza dell’innesto con le due pareti 


perpendicolari al piano dell’arco. 
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Figura 38: Posizione 
5 
dell'aro trionfale nella 
pianta della chiesa. 


VAASCICLA 
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SL 


x TRAGSI 
| en 


ù) 
Figura 39: Modello FEM utilizzato per l'analisi statica non lineare. 


10,6 


46 d7 


In figura 40 viene rappresentata la mappa delle 
lesioni valutate in corrispondenza del carico 
massimo (analisi statica non lineare), che 
suggerisce la posizione delle cerniere che 


caratterizzano il meccanismo di collasso. ma 


Figura 40: Crack pattern in corrispondenza del carico 


Risultati e conclusioni 


Sulla base del confronto integrato fra analisi del danno, meccanismi di collasso proposti dalle linee 


guida e modellazione agli elementi finiti, si propongono le seguenti osservazioni conclusive: 


e Imeccanismi di crisi suddivisi per macro-elementi appaiono idonei in termini generali come 
approccio semplificativo allo studio della risposta al sisma di strutture monumentali, ma sono 
in assenza di sensibili asimmetrie. 

e L’analisi modale con il FEM in campo elastico appare rispondere ben rispetto ai meccanismi 
di crisi che si sono poi effettivamente attivati a conferma del fatto che detta procedura è 
consigliabile per livelli di conoscenza complessivi ma non sempre di dettaglio; in particolare, 
il modello coglie bene gli effetti indotti dalla navata laterale e il meccanismo di danno 


evidenziato dal campanile. 
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2.3.1.5 Basilica di Santa Maria all’ Impruneta 


STATO Italia 
CITTÀ Impruneta (FI) 
ANNO DI 
COSTRUZIONE XIV secolo 
https://www.researchgate.net/publication/ 
da 235758858 
Tabella 11: Identificazione chiesa. Figura 42:Vista dall'alto, Google Maps. 


Figura 41:Basilia di Santa Maria all’Impruneta. 
Cenni storici e descrizione della struttura 


La disposizione attuale della Basilica risale al XIV secolo 
come ampliamento di una piccola chiesa romana consacrata } 


nel 1060, su un santuario etrusco-romano, dovuto alla 


crescente fama di un'immagine miracolosa della Madonna. 


L’abside poligonale e le due cappellette ai lati del transetto 


sono stati costruiti nel XVI secolo. Nonostante molte 


alterazioni, distruzioni e restauri nel seguito dei secoli, 


l’impianto trecentesco è quello che vediamo ancora oggi. 


40m 


Caratteristico della basilica è il pronao costruito nel 1634, che 


è composto da cinque ampie arcate, sormontate da finestre 


rettangolari che ricoprono quasi totalmente la vecchia 


facciata. 


La vista in pianta (figura 43) mostra la tipica pianta a navata 


unica con un’abside poligonale, separato da un arco trionfale. 


6m 


Le dimensioni principali sono 39,90 m di lunghezza, 14,50m | "È 
Figura 43: Planimetria della chiesa. 
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di larghezza massima e un’altezza del muro di circa 15 m. Lo spessore delle pareti varia da 0,70 m 
(navata) e 0,80 m (abside) e la copertura è una struttura in legno. Differenti tipologie di muratura si 
possono notare, differenti sia per materiale (pietra, mattoni ecc.) che per forma (“opus incertum”, 
“opus mixtum” ecc.), anche se la costruzione è prevalentemente in muratura di arenaria irregolare 
(pietra locale) con giunti di malta spessa. Siccome la chiesa non presenta un particolare stato di danno, 
non è stata fatta una completa caratterizzazione dei materiali, ma sono state prese informazioni già 


disponibili di materiali simili. 
Modellazione FEM e analisi svolte 


La vulnerabilità sismica della chiesa è stata valutata con un approccio combinato tra due metodi di 
analisi. Il primo è un’analisi globale tramite la costruzione di un modello agli elementi finiti (figura 
44), che tiene conto del comportamento non lineare della muratura utilizzando specifici legami 
costitutivi, mentre il secondo è un’analisi limite semplificata su alcuni macro-elementi, scelti in 


riferimento ai meccanismi attivati su edifici simili durante i terremoti passati. 


Software Ansys v.11.0 

Tipologia di Elemento *Pareti: elemento solido isoparametrico con 8 nodi + Solid65 

*Volta: elemento bidimensionale isoparametrico con 4 nodi + Shell143 
*Capriata a traliccio: elemento unidimensionale isoparametrico con 2 
nodi >Beam44 

Dimensione *27779 nodi; 

‘76895 elementi 3D; 

*1751 elementi 2D; 

*547 elementi 1D. 

Utilizzo *Valutazione comportamento strutturale statico e vulnerabilità sismica; 


Tabella 12: Descrizione modello FEM. 


Figura 44: Modello agli elementi finiti della Basilia di Santa Maria all’Impruneta 
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Un'indagine preliminare in situ è stata fatta per ottenere informazioni base sulla geometria, dettagli 
strutturali e irregolarità. Per rappresentare il comportamento non lineare è stato utilizzato il legame 
costitutivo perfettamente plastico Drucker-Prager e il criterio William e Warnke,(Fig. 45) per 
rappresentare la superficie di rottura. Il modello è considerato incastrato alla base. 

Il comportamento statico strutturale è stato analizzato sotto dei carichi verticali, derivanti dal peso 
proprio e dai carichi della copertura (figura 46) ed è stata eseguita anche un’analisi modale (figura 


47). 


/ 


Figura 45: Dominio Druker-Prager e 
William-Warnke 


-1.274 
-.837706 
-.401515 

.034676 
+470867 


+907057 


Figura 47: Tensione verticale oz (N/mm2), vista globale Figura 46: Primo modo di vibrare della chiesa. 
Mode no. Period (s) X-direction (transversal) Y-direction (longitudinal) Z-direction (vertical) 
Meg (*) EM (%) Mar(%) EM (%) Meg (%) EMyr(%) 

1 0.585 203 10.58 14.11 23.95 0.00 024 
2 0.480 292 13.50 001 23.96 0.00 024 
3 0.404 002 13.53 072 2468 0.00 025 
4 0.385 8.62 22.15 128 25.96 0.00 025 
5 0334 0.93 27.57 1.13 27.68 0.01 036 
6 0298 019 27.76 18.38 46.06 0.26 062 
7 0297 202 27.78 531 51.37 2.86 348 
8 0289 073 28.53 0.67 55.22 0.00 843 
9 0261 610 34.64 0.00 55.22 0.28 871 
10 0242 0.29 35.11 5.66 6287 0.12 8.88 


Tabella 13: Masse effettive modali per la direzione trasversale, longitudinale e verticale. 
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Figura 48: Pushover, cracking pattern. Carico uniforme in direzione y (sinistra), carico uniforme in direzione x (destra). 


Per valutare il comportamento sotto l’azione sismica è stata effettuata un’analisi pushover, dove si 


può valutare lo stato di danno che la struttura può sviluppare sotto l’azione sismica (figure 48). 


I risultati ottenuti con la modellazione agli elementi finiti possono essere confrontati con quelli 


derivanti dall’analisi semplificata di meccanismi di collasso di macro-elementi costituenti la struttura. 


Si individuano diversi meccanismi possibili, per lo più fuori dal piano, grazie all’analisi globale 


precedentemente effettuata. (figura 49 per esempio). 


Figura 49: Meccanismo 2 (locale). Ribaltamento del pronao. 


Mechanism 1 (pronaos) 
Mechanism 2 (pronaos) 
Mechanism 3 (pronaos) 
Mechanism 4 (lateral wall) 
Mechanism 5 (lateral wall) 
Mechanism 6 (lateral wall) 
Mechanism 7 (nave) 
Mechanism 8 (aspe) 


Collapse 
multiplier (%o) 


0.19 
0.25 
0.20 
0.06 
0.64 
0.14 
0.07 
0.10 


Seismic spectral 
acceleration (x) (g) 


0.22 


Seismic 
demand (g) 


0.38 
0.39 
0.33 
0.39 
0.39 
0.29 
0.36 
0.36 


Figura 14: Moltiplicatore di collasso. 
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Risultati e conclusioni 


Lo scopo di questo studio è quello di dimostrare che le analisi numeriche avanzate possono offrire 
importanti informazioni sulla comprensione dell’effettivo comportamento strutturale degli edifici 
storici in muratura. Si ritiene che la metodologia usata e le conclusioni di questo caso studio siano 


applicabili a un’ampia varietà delle Basiliche esistenti nel territorio. 


2.3.2 Casi studio FEM con Model Updating 
2.3.2.1 Chiesa gotica di Fossanova 


STATO Italia 
CITTÀ Priverno (LT) ; 
ANNO DI COSTRUZIONE 1187 d.C. \ Pol 
iScotteria Del\BOrgo N % Da 
de 
FONTE sc... da 


Tabella 15: Identificazione chiesa. 


Figura 51: Chiesa gotica di Fossanova 


Cenni storici e descrizione della struttura 


L'Abbazia di Fossanova nacque verso la fine del ventesimo secolo dalla trasformazione di un 
preesistente monastero benedettino, forse risalente al sesto secolo, di cui rimane una minima traccia 
al di sopra del rosone della chiesa. Terminata nel 1205 e consacrata nel 1208, è considerata il primo 
monumento nazionale (dichiarata nel 1874) dell’architettura primitiva a tutto sesto introdotta in Italia 
dai Cistercensi Francesi. Il complesso architettonico fu probabilmente progettato e realizzato dai 


francesi e dai locali monaci. 
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Dalla figura 52 si nota che è composta da tre navate, con sette 
campate, transetti e abside rettangolare. Tra la navata principale 
e il transetto si trova un’alta torretta con il campanile. Le 
dimensioni principali sono 69,9 m di lunghezza, 20,1 m di 
altezza e 23,2 metri di larghezza. 

La navata, i corridoi, il transetto e l'abside sono coperte da volte 
a tutto sesto. 


Nel corso dei secoli il complesso ha subito alcune piccole 


modifiche estetiche nel prospetto principale e in copertura. 
Sono state svolte prove in situ (es. endoscopia), e di laboratorio 


(es. prova di compressione e analisi microscopica sulla malta) 


per la caratterizzazione meccanica dei principali elementi e il 
materiale base usato per la chiesa è un calcare sedimentario dci 


molto compatto. 


Prova di vibrazione ambientale 


Le caratteristiche dinamiche dell’intera chiesa sono state stimate da prove di vibrazione ambientale, 
con strumenti in diversi punti della facciata, delle volte e delle tre navate. Le analisi sono state 
eseguite nelle due direzioni principali ortogonali, considerando tanti punti preselezionati lungo la 
lunghezza delle colonne e dei contrafforti. Inoltre, sono state effettuate delle misure nella direzione 


verticale sugli archi e le volte del piano del tetto. 


Relat. amplitude 
Relat. amplitude 


0 25 


(a) frequency (Hz) (b) frequency (Hz) 


Figura 53: Grafico dello spettro amplificato di Fourier: a) direzione longitudinale; b) direzione trasversale. 


I risultati ottenuti mostrano le frequenze fondamentali della chiesa nelle direzioni principali 
orizzontali. Come evidenziato in figura 53, lungo la direzione trasversale il primo picco ha una 
frequenza di 3,8 Hz (primo modo in direzione trasversale), mentre lungo la direzione longitudinale il 


primo picco ha una frequenza di 4,6 Hz (primo modo in direzione longitudinale). Il primo modo 
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torsionale, dai risultati mostrati in entrambe le direzioni, si trova tra 6,6 e 6,8 Hz. Il range di frequenza 


tra 2,4 e 2,8 non deve essere considerato, poiché queste frequenze sono dovute a vibrazioni esterne. 


Modellazione FEM e analisi svolte 

Le analisi svolte hanno l’obiettivo di andare a studiare il comportamento sismico della chiesa in 
esame. Dopo la caratterizzazione meccanica dei materiali, sono state fatte delle prove di vibrazione 
ambientale ed è stato costruito anche un modello agli elementi finiti. Il modello FEM, mostrato in 
figura 54, è stato raffinato e calibrato sui risultati provenienti dalla prova di vibrazione ambientale, in 
modo da riprodurre fedelmente il comportamento dinamico della struttura. 

Il modello FEM riproduce la geometria della chiesa, in cui i vari elementi che costituiscono la chiesa 


sono stati modellati tramite elementi solid, shell e beam. 


Software Abaqus 
Tipologia di Elemento |*Solid 
*Shell 
*eBeam 
Dimensione Struttura in muratura Struttura in legno (copertura) 
*59219 nodi; «4759 nodi; 
«217571 elementi solid; *4447 elementi solid; 
*4976 elementi shell. *1987 elementi shell. 
Utilizzo *Valutazione comportamento strutturale e vulnerabilità 
sismica; 
*Confronto con altra tipologia di modellazione (shaking- 
table test su modello a scala 1:5,5). 
Calibratura Prova di vibrazione ambientale 


Tabella 16: Descrizione modello FEM 


Figura 54: Modellazione agli elementi finiti della chiesa. 
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Model Updating 


x 


Per prima cosa è stata effettuata un’analisi modale per determinare le frequenze naturali 


corrispondenti dei principali modi di vibrare. Sono stati sviluppati due modelli. 


e FEMI: calibrato rispetto alle frequenze ottenute sperimentalmente, cambiando il valore del 
modulo di Young dei pilastri delle volte e delle pareti (riferito a un equivalente modulo di 
Young per gli edifici in muratura). 

e FEM2: evoluzione del modello FEMI. Non vengono considerate le strutture in legno della 


copertura, per valutare la loro influenza. 


Dalla tabella 17 si nota che le frequenze del modello FEMI sono molto vicine alle frequenze date 
dalle vibrazioni ambientali. Questo significa che il modello proposto è in grado di riprodurre la 


risposta dinamica della chiesa. 


Detected natural frequency value (Hz) 


Direction Ambient vibration Finite element analysis (FEM) (1)* FEM (2)** 
Transverse 3.8 3.77 3.65 
Longitudinal 4.6 4.63 4.30 
Torsional 6.8 6.45 6.15 


*With roofing wooden elements. 
**Without roofing wooden elements. 


Tabella 17: Confronto delle frequenze naturali. 


Figura 55: Modi di vibrare: a) trasversale, b) longitudinale, c) torsionale. 


Una volta ottenuto il modello sono state svolte delle analisi non lineari che hanno consentito 


l’identificazione delle parti più vulnerabili della chiesa (figura 56). 
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8, 833 
(ave. Crit.: 758) 


Ss, S33 
(Ave. Crit.: 75%) 


241e-01 +3.521e+01 
+3 .,000e-01 +9.000e-01 
+1.083e-01 +8 .250e-01 
-8.333e-02 +7.500e-01 
-2.750e-01 +6.750e-01 
-4.667e-01 +6.000e-01 
-6.583e-01 +5.250e-01 
-8.500e-01 +4,500e-01 
-1.042e+00 +3.750e-01 
-1.233e+00 +3.000e-01 
-1.425e+00 +2 sneni 
-1.617e+00 +1 e-D1 
-1.808e+00 +7.500e-02 
-2.000e+00 +0.000e+00 

(a) (b) 


Figura 56: a) analisi statica con i pesi propri, b) analisi dinamica spettrale, con la relativa distribuzione di sforzo. 


Prendiamo per esempio una sezione delle tre navate (figura 57), dove si raggiungono alti valori di 


trazione. 


progressive 
cracking 


® 0.046g 


Le cerniere presenti sono. 
: . . 0.052g 
e  duenegli archi centrali (h=20 m); 


e quattro negli archi adiacenti (h=8 m); sit 


0.066g 


() 
(©) 
0 
(h=0 m). Ò 0.151g 


Figura 57: Parti chiave soggette al più grande sforzo. 


e quattro alla base dei pilastri e contrafforti 


Risultati e conclusioni 


La valutazione della risposta sismica della chiesa di Fossanova è valutata per mezzo del modello FEM 
e delle analisi limite. Infatti, il modello FEM consente di determinare la riposta globale dell’intera 
struttura, prendendo in considerazione l’interazione tra le diverse parti che la costituiscono, e la 
localizzazione delle parti in cui avviene la rottura. Con le analisi limite possiamo avere una stima 
dell’intensità sismica che produce la rottura della struttura secondo i cinematismi previsti. 

Il modello FEMI riproduce fedelmente il comportamento della struttura, grazie alla calibratura 
rispetto ai dati provenienti dalle prove di vibrazione ambientale. Questo passaggio è stato 
fondamentale, poiché ha mostrato come la presenza della struttura in legno di copertura influenza 


effettivamente il comportamento. 
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2.3.2.2 Campanile della Cattedrale Metropolitana di Fermo 


STATO Italia 
CITTÀ Fermo (FM) 
ANNO DI 
COSTRUZIONE 1227/46, 
https://www.mdpi.com/1424- 
FONTE 8220/20/11/3315 


Tabella 18: Identificazione chiesa. 


Figura 59: Cattedrale Metropolitana di Fermo 
Cenni storici e descrizione della struttura del campanile 


Il soggetto dello studio è la Cattedrale Metropolitana di 
Fermo, dedicata a Santa Maria Assunta in Cielo, dove si 
pone l’attenzione sul comportamento dinamico del 
campanile. 
Le prime notizie della sua esistenza risalgono 
probabilmente al secondo secolo e nel corso degli anni ha 
subito modifiche fino ad arrivare alla corrente 
configurazione (figura 60). Nel 1176 la chiesa esistente fu 
distrutta da un incendio causato da un assedio alla città e 
nel 1227 la chiesa venne ricostruita in stile Romano- 
Gotico. Successivamente la chiesa venne demolita e 
ricostruita tra il 1781 e il 1789, preservando la facciata e il 
campanile. La particolare asimmetria della facciata è 
dovuta alla costruzione del campanile in corrispondenza 
della navata della chiesa. Questo comporta che il portone 
e il rosone non si trovano esattamente al di sotto della parte 


più alta della facciata. Nel 1898 vennero predisposti altri 
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Figura 60: Planimetria ed evoluzione della chiesa 


interventi di restauro riguardanti il campanile. Dopo il terremoto del 2016 la chiesa fu danneggiata 


ed è stata chiusa per più di un anno per la restaurazione. 


Il campanile, mostrato in figura 61, è collegato a sud con la facciata e ad est con la chiesa ed è alto 
48,1 m con una sezione quadrata 10x9,5 m. Le pareti sono alte 24,3 m con uno spessore di 2,7 m, 
dove all’interno è presente un vano scala. Il quarto e il quindi piano hanno muri di spessore ridotto 
(1,5 m) ed è presente una finestra su ogni lato, costituita da due aperture ad arco che poggiano su due 
colonne. L’ultimo piano è a forma ottagonale, dove su ciascun lato, lungo circa 3,1 m, è presente una 
finestra circolare e le pareti hanno uno spessore di 0,5 m. Gli ultimi tre livelli sono collegati da una 
scala a chiocciola situata nell’angolo nord-est. I muri sono costituiti da mattoni pieni e pietre e i solai 


sono costituiti da volte in muratura. 


North front West front Sud front East front Section X-X 


A. x XxX XxX XxX XxX XxX XxX x 


Ground floor 1 floor 2"4 floor 3"4 floor 4° floor 5°" floor 6°" floor 


Figura 61: Caratteristiche geometriche del campanile. 


Monitoraggio della salute strutturale (SHM) 


Il monitoraggio SHM ci fornisce il comportamento dinamico reale della struttura. Questo è 
fondamentale per andare a calibrare il modello numerico visto che restituisce dei valori diversi da 
quelli reali a causa di diversi aspetti, come per esempio approssimazioni geometriche, connessioni 
tra gli elementi, assunzioni riguardo i parametri dei materiali. Per superare questi problemi vengono 
usati i risultati OMA provenienti dal monitoraggio delle vibrazioni della struttura prodotte da un input 
sconosciuto (vibrazioni ambientali) sotto determinate condizioni operative. Il posizionamento dei 
sensori viene mostrato in figura 62, dove sono stati fissati con una resina in modo da connetterli 


correttamente alla struttura. 
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4. 
Figura 62: Posizionamento sensori. 
L’identificazione delle frequenze è stata SSI Method 
0 
implementata usando il metodo SSI (Stochastic fi = 1654 Hz 
20 | 
Subspace Identification), sviluppato nel dominio af” 1800 Hz 
RIOT : N 40 fa = 4236 Hz 
del tempo. La figura 63 mostra il diagramma di È - L 
Ln. : . . o. Ta Il | f,=4955Hz fa = 9.767 Hz 
stabilizzazione usato per separare i modi fisici e è Il | | 
5 II I | 
[m) a Ì 
simulati, con il quale è possibile vedere  ©®| , È NA JM À SRI sella, # 
pit dd VANI. LI Mi Mia vi . 
graficamente le frequenze naturali. 200 | - I 1 Ni 
Il LI I 
120 Il Luo Il 
0 2 4 6 8 10 12 
f 182] 
Figura 63: Grafico di stabilizzazione 
Modellazione FEM 


In questo studio sono stati costruiti due modelli in due software diversi in modo da poter svolgere due 


differenti modalità di aggiornamento e confrontare i risultati ottenuti. 


i MODELLOFEM_C_ 
Software *Midas FEA 
*Code Aster (per l'updating tramite l'algoritmo 
genetico) 
Tipologia di Elemento | Elemento solido tetraedrico con 4 nodi 
Utilizzo Identificazione dinamica 
Calibratura Monitoraggio SHM 


Tabella 19: Descrizione modello FEM 
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Figura 64: Modello agli elementi finiti: (a) in prospettiva, (b) sezione. 


È stato costruito un modello dell’intera struttura poiché la reale interazione tra il campanile e la chiesa 
era sconosciuto. In un primo modello (step 0) sono state assegnate le caratteristiche meccaniche del 
materiale presente considerando il limite inferiore. Inoltre, i parametri elastici assegnati al modello 
sono stati corretti con dei coefficienti presenti nella norma (Tabella.C8.5II presente nella Circolare 
applicativa alle NTC18) a causa di alcune limitazioni come una mera ispezione visiva e la mancanza 
di informazioni sul materiale utilizzato e il suo stato conservativo. Infine, la base è stata considerata 
incastrata a terra con l’inserimento di un vincolo e sono stati assegnati i carichi dai dai solai e dalle 
campane. Nel modello di partenza il modulo di Young è considerato uniforme e uguale a 1800 Mpa. 
Prima di procedere con l’updating del modello, sono state confrontate le frequenze provenienti dal 
modello sperimentale (EM) e dal modello numerico (NM). Questa fare viene identificata come step 


0, e la tabella 20 mostra la grande differenza tra i due modelli, in termini di frequenza. 


Step 0 


Mode fem [Hz] bui Afp (%) 
1 1.654 1.337 19.16 
2 1.800 1.671 7.17 
3 4.236 3.876 8.49 
4 4.955 4.182 15.60 
5 9.767 7.404 24.19 


Tabella 20: Confronto tra le frequenze allo step 0 del modello sperimentale 
(EM) e del modello numerico (NM). 


Model Updating 


L’aggiornamento del modello è stato svolto in due modalità. 


> updating manuale modificando i parametri elastici dei materiali, data la disponibilità di un 


rilievo geometrico ad alta precisione. Un’ispezione visiva ha mostrato che la pietra è 
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concentrata nella prima metà della struttura. Così è stato riprodotto un primo modello, 
dividendo il campanile in sezioni, nelle quali il modulo elastico diminuiva dal basso verso 
l’alto. In totale sono stati assegnati 4 moduli elastici (step 1) aumentando 


significativamente tutte le frequenze eccetto la terza (tabella 21). 


Step 1 


Mode fem [Hz] tag Eb] Afi (%) 
1 1.654 1.500 9.31 
2 1.800 1.903 —5.72 
3 4.236 3.797 10.36 
4 4.955 4.969 —0.28 
5 9.767 8.516 12.81 


Tabella 21: Confronto tra le frequenze allo step I di EM e NM. 


Successivamente il modulo elastico venne nuovamente modificato lungo lo spessore delle 
pareti, considerando che la pietra Istriana era presente solamente nel paramento esterno. 
Inoltre, siccome la porzione della torre nel lato nord appare costituito da materiale 
disomogeneo, il modulo elastico venne ridotto di molto. Queste assunzioni portarono alla 
formazione di un secondo modello (step 2), nel quale sono stati assegnati 23 diversi 
moduli elastici. Dalla figura 65 e dalla tabella 22 si nota che l’errore è stato ridotto 


notevolmente e che si è ottenuta una buona corrispondenza in termini di frequenza. 


Step 2 o 
Mode fem [Hz] tons [x] Afp (%) 

1 1.654 1.610 2.66 

2 1.800 1.871 —3.94 

3 4.236 4.272 —0.85 

4 4.955 5.036 —1.63 

5 9.767 8.504 12.93 


Tabella 22: Confronto tra le frequenze allo step 2 di EM e NM. 


Mode 1 Moae Z Mode 3 Mode 4 Mode 5 
NM fir 1.610 Hz fa= 1.871 Hz far 4.272 Hz fa” 5.036 Hz fs= 8.510 Hz 
Trans. Y Trans. X Tors. Z Flex. Y Tors. Z + Flex. Y 
ì i la gp” » Ì | 
M | M [1 A Ti Mod V ] 
EM fy= 1.654 Hz fs= 1.800 Hz fy= 4.236 Hz fa 4.955 Hz fe 9.767 Hz 
Trans. Y Trans. X Tors. Z Flex. Y Tors. Z + Flex. Y 


Figura 65: Confronto tra le frequenze allo step 2 di EM e NM. 
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Un ulteriore confronto si può valutare con il MAC (modal assurance criterion) e 1 risultati 
in figura 66 mostrano che la CrossMAC tra EM e NM non fornisce valori molto 


soddisfacenti. 


Cross MAC - Step O Cross AC - Step 1 


7 ae > VA? 


CrossM AC . Step 2 


sets. 
to) è fa 


Figura 66: Cross-Modal Assurance Criterion (CrossMAC) tra EM e NM. 


> updating tramite l’utilizzo di un algoritmo genetico. Questo algoritmo considera un range 
di valori possibili dei parametri considerati, dal valore minimo al valore massimo, con i 
quali svolge le diverse analisi andando a trovare le soluzioni migliori. Anche qui sono stati 
sviluppati due modelli nel software Code Aster, dove le proprietà dei materiali sono state 
assegnate sempre a elementi solidi tetraedrici e gli step per l’aggiornamento sono gli 
stessi. 
In primo luogo, è stata svolta un’analisi modale allo step 0, dove i risultati sono 
esattamente uguali a quelli del modello NM del codice Midas FEA. 
Nelle tabelle 24 e 26 vengono mostrati i risultati degli step 1 e 2, dove vengono inseriti 
valori del modulo elastico diversi dalla modalità manuale, poiché l’algoritmo cerca di 
trovare i valori che meglio ottengono la convergenza dei valori delle cinque frequenze 


(Tabelle 23 e 25). 


45 


Masuns Tadtntio 
secca o neo ace cea 
«DRM EZZZAZZIEZEO 


234%» vo 


(b) 


Figura 67: Modelli in base al materiale assunto. a) quattro materiali 
b) ventitre materiali. 


Number of Material E [MPa] Manual Calibration E [MPa] GA Calibration 


1 2700 1000 
2 3200 5090 
3 2800 2230 
2 2700 4540 


Tabella 23: Differenti moduli elastici E [MPa] tra la modalità manuale e automatica allo step 1. 


Mode fem [Hz] iù Aff} (%) 
1 1654 1.645 0.62 
2 1.800 1.962 8.95 
3 4236 4848 1445 
4 4955 5.134 3.61 
5 9.767 8.567 12.29 


Tabella 24: Confronto tra le frequenze di EM e NM allo step 1 con Code Aster. 


Number of Material E [MPa] Manual Calibration E [MPa] G.A. Calibration 


1 1200 3797 
2 2700 3425 
9 1800 4005 
4 1250 1554 
5 1800 1283 
6 1800 1592 
7 4400 4473 
8 3000 1258 
9 2000 2917 
10 3800 1092 
11 2500 1970 
12 1800 2322 
13 800 4096 
14 1500 3828 
15 1200 4137 
16 1800 3588 
17 1800 2865 
18 3200 3550 
19 4400 3249 
20 2800 2251 
21 1900 1309 
22 1800 3295 
23 2400 871 


Tabella 25: Differenti moduli elastici E [MPa] tra la modalità manuale e automatica allo step 2. 
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Step 2 


ly 

Mode feM [Hz] tum [Hz] Af (%) 
1 1.654 1.529 6.97 
2 1.800 2.000 1.95 
3 4.236 3.715 23.36 
4 4.955 5.095 0.76 
5 9.767 9.686 13.05 


Tabella 26: Confronto tra le frequenze di EM e NM allo step 2 con Code Aster. 


In figura 68 sono riportati i risultati della CrossMAC. 


CrossMAC with Code Aster® - Step 1 


CrossMAC with Code Aster© - Step 2 


Figura 68: Cross-Modal Assurance Criterion (CrossMAC) tra EM e NM. A sinistra step 1 e a destra step 
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Risultati e conclusioni 


La ricerca presentata in questo articolo è partita dall’osservazione diretta delle caratteristiche 
dinamiche del campanile della cattedrale metropolitana di Fermo tramite i risultati dati da OMA con 
l’obiettivo di confrontare due diverse procedure di calibrazione: una manuale e una automatica. In 
entrambe le procedure sono state perturbate solo le matrici di rigidità, andando a modificare il modulo 
elastico dei materiali presenti. Analizzando 1 risultati derivanti dalla calibrazione manuale si è visto 
che esiste una corrispondenza valida tra le frequenze numeriche e sperimentali, poiché è possibile 
avere un controllo maggiore sulla distribuzione del modulo di Young, anche se la matrice CrossM AC 
non ha dato risultati soddisfacenti. Migliori risultati sono stati ottenuti con la calibrazione automatica 
nella matrice CrossMAC, mentre le frequenze hanno una deviazione maggiore. Questo può derivare 
da due cose. La prima è che l’aggiornamento con l’algoritmo è stato fatto non solo considerando i 
valori di frequenza, ma anche esaminando il CrossMAC, e la seconda cosa è la possibile influenza 
dell’interazione suolo-struttura delle pareti delle parti che interfacciano con il campanile che non sono 


state direttamente strumentate e che forse potrebbero mostrare una non perfetta interazione. 
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2.3.2.3 Chiesa di San Torcato 


STATO Portogallo 
CITTÀ San Torcato (Guimaraes) 
ANNO DI 
COSTRUZIONE 187146, 
https://\www.researchgate.net/publication/22 
FONTE 6232385 


Tabella 27: Identificazione chiesa. 


Figura 70: Chiesa di San Torcato. 


Descrizione della struttura 


La chiesa di San Torcato si trova nel paese di San 
Torcato, vicino la città di Guimaraes, in 
Portogallo. La sua costruzione iniziò nel 1871 ed 
è stata completata solo recentemente. Dalla 
forma della pianta architettonica (figura 71) si 
nota lo schema a croce classico di una basilica 
con la navata principale coperta da una volta a 
botte, con il transetto sovrastato al centro da una 
cupola e l’abside. Nella facciata sono presenti 


due torri e tra queste e il transetto venne aggiunto 
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Figura 71: Planimetria della chiesa. 


un piano di altezza ridotta in entrambi i lati della navata, dove sono presenti dei locali dai quali è 


possibile accedere alle torri. La navata principale ha come dimensioni principali 25 m di lunghezza e 


15 m di larghezza, con i muri di 1,30 m di spessore, e la volta a botte poggia su archi in pietra, che 


scaricano su colonne in pietra, supportate da dei contrafforti adiacenti ai muri della chiesa. La cupola 
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principale, all’interno del transetto, è sopportata da un doppio arco in muratura. Le torri presentano 
una sezione rettangolare di 7,50 x 6,20 m, misurata dalla superficie esterna, con le pareti spesse 1,45 
m un’altezza totale di 32,70 m. A causa della lunga durata della costruzione, le tecniche strutturali e 
i materiali variano all’interno dell’intera struttura. Il materiale principale e originario è la pietra 
granitica, posata con giunti a secco, ma è presente anche qualche parete in cemento armato. Al 
momento la chiesa presenta dei problemi strutturali dovuti al cedimento del terreno nella facciata 
principale. La continuità delle lesioni all’interno della chiesa e l’inclinazione delle torri indica un 
assestamento dovuto all’elevato libello di stress delle torri e al soffice strato di riempimento del 
terreno. Per stabilizzare il danno è previsto un intervento strutturale ed è già stato predisposto un 


monitoraggio per assistere l'intervento. 


In questo studio sono stati effettuati test di identificazione modale per andare a calibrare il modello 
agli elementi finiti, utilizzato per le analisi della sicurezza. Il lavoro comprende la campagna 
sperimentale per la stima dei parametri modali (frequenze naturali, forme modali e coefficienti di 


smorzamento) e procedure per calibrare e aggiornare accuratamente il modello FE. 
Identificazione dinamica 


Le analisi sperimentali modali sono state svolte tramite l’utilizzo di accelerometri e le misure sono 
state prese in 35 punti divisi in 9 prove. Gli accelerometri di riferimento sono stati posizionati in cima 
alle torri a causa della loro ampiezza elevata e del contributo modale in ogni modo. Le torri e la 
facciata sono state misurate accuratamente visti i loro seri danni, così il numero e la posizione dei 
punti di misura sono prevalentemente presenti in questa area (figura 72). Come tipologia di 


eccitazione è stata scelta la vibrazione ambientale. 


Figura 72: Posizionamento sensori. 
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Con la tecnica FDD (Frequency Domain Decomposition) ogni modo è stato stimato come la 
decomposizione delle densità spettrali della risposta del sistema in diversi sistemi ad un singolo grado 
di libertà (SDOF) e inoltre i risultati sono stati confrontati con altre tecniche. I risultati dei primi 


quattro modi di vibrare sono visibili in figura 73 e tabella 28. 


DICE 9F/Hz] Frequency Domain Deccmpcation - Peak Ackng 
Average of Ue Nov mozioni Simpad ee Voluso at rpricolira EA 
Spertrel Density Melrices ot al Test Setups reqrey = 2 


D Estates Moce 
+ Salociod Mode 
== Node Marher 
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ND Lire 31 
SVD Une 3. 
SVD Line 33 
SND Lino #4 
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[Dmera nni Ratiti earn i Faina e 5pa0S MORI Sato Space Dimeetine 
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“ Indicatore 
an MIE Curi 
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i) A SKD Subopore 
PA ® Sil 
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40 — — Mode Marker 
sl, Unes 
b.] YD Lire 81 
ii SAD line #2 
— S\D Line 33 
SVD Line #6 
Frequency ha] 
Figura 73: Metodi di stima: a) metodo FFD, b) metodo SSI. 
FDD EFDD CFDD SSI 
Frequency — Error Frequency Error Frequency Error [%] Frequency 
Hz % Hz % Hz Hz 
Mode 1 2.15 04 2.14 0.1 2.13 0.3 2.14 
Mode 2 2.64 04 2.62 0.1 2.62 0.1 2.63 
Mode 3 2.89 13 2.89 1.1 2.86 0.3 2.85 
Mode 4 2.97 14 2.94 0.3 2.93 0.1 2.93 


Tabella 28: Confronto tra i metodi del dominio di frequenza con il metodo SSI. 


Modellazione FEM 


Sono stati modellati solamente la facciata, la navata e le torri (figura 74). Il suolo è definito elastico 
e gli elementi in muratura sono considerati omogenei in tutte le parti della struttura. 


Wi hh 


Figura 74: Modello agli elementi finiti della chiesa. 
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Software iDiana Release 9.3 


Tipologia di Elemento | Struttura Interazione con terreno 
«Elemento solido Elemento di intercaccia 
quadratico con 20 nodi> |quadrilatero con 16 nodi> 
CHX60; CQ48. 


«Elemento wedge 
quadratico con 15 nodi 


CTP45. 
Dimensione *3044 elementi 

*3153 nodi 
Utilizzo Identificazione modale 
Calibratura Monitoraggio dinamico 


Tabella 29: Descrizione modello FEM. 


Le seguenti ipotesi e assunzioni sono usate nel modello: 


l’interazione suolo-struttura è definita con elementi di interfaccia che consentono la 
definizione delle proprietà di rigidezza verticale e orizzontale per il suolo. Per le assunzioni 
delle proprietà, sono state usati i risultati di investigazioni precedenti sul suolo; 

Peso della muratura di 25 KN/m}; 

Il coefficiente di Poisson fissato a 0,2; 

In tutta la struttura è usato il materiale muratura omogeneo; 

A causa del modello incompleto, la parte mancante della struttura è stata simulata con 
l’introduzione degli elementi di interfaccia all’intersezione con i transetti, per simulare la 


rigidezza delle parti mancanti. 


Model Updating 


Vengono aggiornati 4 parametri (Em, Es, En, Ets) e 1 modi sperimentali e numerici vengono confrontati 


tramite il Modal Assurance Criterion (MAC), che è il metodo più comunemente usato per investigare 


sulla coerenza dei due differenti modelli. Dipende dalla correlazione dei vettori modali, varia tra 0 e 


1 ed è definito da 


Liza gÎ gi l° 


MACen= ST 7AaSSrTTAG 
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dove g° e p" sono i vettori dei modi di vibrare, sperimentali e numerici rispettivamente. 
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MAC & Frequency Comparison 
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Figura 75: Confronto tra i risultati numeri e sperimentali prima dell'updating. 


Inizialmente il contributo di ogni parametro è stato studiato indipendentemente dagli altri e in figura 
75 possiamo vedere il risultato del primo step per l'aggiornamento del modello in termini di confronto 
tra le forme modali, in base alla parte della struttura considerata, e il confronto tra i valori MAC e 
della frequenza. 

Nel secondo step sono stati modificati i parametri considerati, utilizzando il metodo Douglas-Reid, 
dove i limiti inferiori e superiori per i valori sono stati definiti secondo il giudizio ingegneristico 


(tabella 30). 


Name Optimum Modes ex [Hz] @re[Hz] @Error[%] MAC 
Em 5.642 GPa 1° (Transversal Y ) 2.14 2.14 0.09 0.92 
Es 0.629 GPa 2°“ (Longitudinal X) 2.63 2.55 2.86 0.86 
En 0.046 GPa 3" (Torsional) 2.85 2.93 2.70 0.83 
Ets 21.591 GPa 4° (Torsional) 2.93 2.94 0.41 0.77 


Tabella 30: Risultati del modello ottimizzato. 
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Experimental Results FE Model 
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Figura 76: Confronto tra i modi di vibrare 


I modi di vibrare mostrano una buona coerenza con i modi sperimentali (figura 76). Nel primo e 
secondo modo (modi traslazionali), nessuna differenza significativa è osservata tra i modi. Il terzo 
modo di vibrare del modello numerico presenta un comportamento simile con quello sperimentale, 
ma se si pone l’attenzione alla torre nord, si può vedere che l’effetto torsionale è più accentuato 
rispetto al modello numerico. Anche il quarto modo mostra un comportamento simile, ma presenta 


un basso valore di MAC a causa dello stato di danno che non è considerato nelle analisi. 
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Figura 77: Risultati final updating: (a) confronto tra MAC e frequenze; (b) evoluzione delle analisi di updating. 
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Risultati e conclusioni 


Lo studio presenta l’analisi di identificazione modale della chiesa di San Torcato, insieme a un’analisi 
di aggiornamento di un modello FE. I dati raccolti sul campo sono stati elaborati e controllati con 
diverse modalità. A causa del basso livello di eccitazione, sono state stimate solamente i primi quattro 
modi di vibrare della struttura. I risultati ottenuti sono di elevata coerenza con 1 risultati preliminari 
in termini di frequenza e forme modali. Un modello FE, precedentemente costruito per l’analisi della 
sicurezza, è stato aggiornato utilizzando procedure di ottimizzazione manuali e robuste, modificando 
le proprietà elastiche omogenee della muratura, del suolo e dei parametri di rigidità degli elementi di 
interfaccia. È stata ottenuta un’elevata coerenza del modello numerico con modifiche manuali e 
algoritmi non lineari, ma per ottenere risultati migliori è opportuno effettuare ulteriori analisi 
considerando il danno. La chiesa continua ad essere osservata da un sistema di monitoraggio dinamico 
continuo, utile sia per valutare l’intervento strutturale programmato, sia per valutare l'evoluzione del 


suo comportamento. 


2.3.2.4 Cattedrale Metropolitana di Santiago 


STATO Cile 


l'agire 


A 


Uateara! 


x . Cc Ò = =" 0 tedrallMetropolitan 
CITTA Santiago = : _ a ‘ de Santiago 
ANNO DI 
COSTRUZI 1746 d.C. 
ONE 


https://www.sciencedirect.com/scie 
FONTE nce/article/abs/pii/S0141029617311 
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Tabella 31: Identificazione chiesa. 


Figura 79: Cattedrale Metropolitana di Santiago 
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Cenni storici e descrizione della struttura 


La Cattedrale Metropolitana di Santiago è un’importante struttura storica, infatti nel 1951 venne 
dichiarata un Monumento Storico Nazionale. La costruzione, iniziata nel 1946, è partita dalla facciata 
occidentale per proseguire verso quella orientale, e nel 1830 terminò la prima fase della costruzione, 
senza ne torri ne cupola. Nel 1858 riprese la costruzione della cappella e di una torre, che non 
corrispondeva a nessuna torre esistente nell’attuale cattedrale. Dopo due terremoti avvenuti nel 1851 
e 1874, la chiesa subì gravi danni e così nel 1897 il vescovo ordinò la ristrutturazione che portò alla 
costruzione delle due torri nella facciata principale, la cupola sopra l’altare, il tetto, e altri 
rimodellamenti architettonici. Nel 1898 fu demolita la torre e i lavori furono completati nel 1906.11 
sistema strutturale è costituito principalmente da pareti in muratura in pietra fino ad un’altezza di 11 
m, poi muratura in laterizio fino ad un’altezza di 17 m. La parte centrale e la navata nord hanno una 
copertura costituita da capriate in legno, mentre nella navata sud ci sono capriate metalliche. Nella 
facciata principale si distinguono due torri campanarie alte 45 m, composte da colonne in acciaio agli 
angoli, coperti da muratura in mattoni. La chiesa si divide in tre navate: quella principale ha un’altezza 
di colmo di 20 m, mentre quelle laterali sono alte 15 m. Nella navata principale (tra gli assi 10 e 11 e 
B e C, figura 80) è presente la cupola principale alta 50 m, dove alla base presenta delle travi in 


acciaio che fanno da cordolo per le pareti sottostanti in muratura. 


110. m APROX, 


è 9 ® 9 o ? © o o 9 © è 
6 ut — È 


ell Bisi 
IU Wy w uw 


(13) Ei 


NORTH LONGITUDINAL have TOWER 7}| 


off s-9-sd_é o_o o_o È 


i I DOME CENTRAL LONGITUDINAL NAVE 
©Chg_ dP P_ VO P_P _Sb_@_G_ __q cu 


SOUTH LONGITUDINAL NAVE TOWER / 3 


mi n 
“ h th ff m ehhh rh sh rh ; 
OLI l i 7 a I 7 T = I 


Ul 
TABERNACLÉ È 
ba all 


SACRISTY A ) 
\ Ì { 
° CHAPEL | 

N AL 


Figura 80: Planimetria della chiesa. 
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Prove di vibrazione ambientale 

Attualmente lo stato del danno presente non compromette la stabilità globale, ma potrebbe interessare 
la stabilità locale di qualche elemento a seguito di un evento sismico. Sono state eseguite delle prove 
di vibrazione ambientale, con le quali andare a calibrare il modello FEM, in modo da definire al 


meglio lo stato del danno e il conseguente comportamento statico strutturale e sismico. 
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Figura 80: Posizionamento dei sensori. 


La figura 17 mostra la posizione degli strumenti di misurazione, dove: 


e SA: base degli archi. 
e ST:inalto sulle pareti. 


e SW: base delle finestre. 


Sono stati usati due differenti metodi: EFDD e SSI. I risultati (modi di vibrare e relative 


frequenze) sono riportati tabella 32. 


Experimental mode EFDD method SSI method 
Frequency (Hz) Frequency (Hz) 
1 1.81 1.74 
2 1.96 1.93 
3 2.03 2.06 
4 2.38 2.38 
5 2.55 2.56 
6 2.76 2.79 
7 2.86 2.86 
8 2.95 2.93 
9 3.18 3.20 


Tabella 32: Confronto delle frequenze ottenute con il metodo EFDD e SSI. 


Modellazione FEM 


Viene costruito un modello FEM (Fig. 81), dove le proprietà meccaniche subiscono due processi di 
ottimizzazione: il primo riguarda anche le condizioni al contorno, date dalle strutture adiacenti, ed è 
basato sulla minimizzazione dell’errore tra le frequenze sperimentali e analitiche, mentre il secondo 
è basato sulla riduzione al minimo della differenza tra frequenza e valori MAC. 


Alla base del modello vi sono le seguenti assunzioni: 
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e Il comportamento della struttura viene considerato come continuo, senza discontinuità 
all’intersezione dei muri. 

e lacupolaele torri campanarie vengono modellate come semplici elementi frame, poiché è la 
struttura rigida ad essere al centro dello studio. 

e Le strutture adiacenti (sacrestia, tabernacolo e cappella) che interagiscono con la struttura, 


vengono modellati come semplici elementi elastici. 


Software *eGID 
*DIANA 

Tipologia di Elemento |Elementi solidi tetraedrici quadratici a 10 nodi: 
*CTE10 

Dimensione «271106 nodi 
*157480 elementi 

Calibratura Prova di vibrazione ambientale 


Tabella 33: Descrizione modello FEM. 


Figura 81: Modello agli elementi finiti della cattedrale 


Model Updating 


La prima fase di ottimizzazione del modello consiste nel confrontare le frequenze ottenute 
sperimentalmente (tabella 34) con le frequenze calcolate con il FEM, basate su dei valori iniziali di 
parametri tabella 36. La corrispondenza tra i modi per l’aggiornamento del modello, è basato sulla 
valutazione qualitativa visiva delle forme modali in figura 82 e sulla quantificazione del MAC (Modal 
Assurance Criterion) tra 1 modi sperimentali e quelli ottenuti dal modello numerico (Tabella 35). Per 
poter essere confrontati, i modi sperimentali richiedono la trasformazione a valori reali (importante 
per la dissipazione di energia). 

I primi quattro modi non vengono considerati poiché sono collegati al movimento delle torri, come i 


modi 9 e 11, poiché dovuti al movimento della cupola. 
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Numerical Frequency Mass participation —Mode description 


mode (Hz) based on FE model 
(3) 

1 1.38 0.48 Bending of south tower. 

2 1.39 0.15 Bending of north tower. 

3 1.41 0.61 Bending of south tower. 

4 1.43 0.43 Bending of north tower. 

5 1.90 27.48 Transverse bending of 
longitudinal walls. 

6 2.02 0.27 Transverse bending of 
south longitudinal walls. 

1) 2.19 1.73 Transverse bending of 
longitudinal walls and 
bending dome. 

8 2.27 0.11 Bending dome. 

9 2.32 155 Transverse bending of 
longitudinal walls and 
bending dome. 

10 2.40 0.53 Transverse bending of 
south longitudinal walls. 

RI 2.54 5.44 Bending dome. 

12 2.56 0.02 Transverse bending of 
longitudinal walls. 

13 2.59 0.63 Transverse bending of 
longitudinal walls. 

14 2.83 2.24 Transverse bending of 
longitudinal walls. 

15 2.97 0.22 Transverse bending of 
longitudinal walls. 


Tabella 34: Frequenze numeriche con i parametri iniziali di calibrazione. 


fex= 2.375 Hz fex3= 2.95 Hz 


fo) : i 7 2 n) i Up Po 
Un ANG RG RG 
fi= 1.898 Hz f.= 2.185 Hz fs= 2.827 Hz fa= 2.965 Hz 
Figura 82: Modi sperimentali e numerici. 
Experimental Mode number Experimental frequency (Hz) Numerical Mode number Numerical frequency (Hz) MAC 
ri. 1.96 5 1.90 0.85 
4 2.38 vd 2.19 0.61 
8 2.95 14 2.83 0.55 
9 3.18 15 2.97 0.69 


Tabella 35: Confronto tra le frequenze ottenute sperimentalmente e numericamente. 
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In tabella 37 si evidenziano i parametri considerati per l’aggiornamento del modello: 
-modulo di Young per la muratura in pietra; 

-modulo di Young per la muratura in mattoni; 

-modulo di Young per la muratura armata; 

-fattore di rigidità laterale dalla Cappella; 

-fattore di rigidità laterale dalla Sacrestia; 


-fattore di rigidità laterale dal Tabernacolo. 


Calibration parameter Initial value 
Young's modulus of stone masonry E.,m (MPa) 1390 
Young's modulus of brick masonry Epm (MPa) 1780 
Young's modulus of the reinforced masonry Erm (MPa) 1780 
Lateral stiffness factor from Chapel 1.00 

Lateral stiffness factor from Sacristy 1.00 

Lateral stiffness factor from Tabernacle 1.00 


Tabella 36: Parametri iniziali. 


Questi parametri di calibrazione hanno un’elevata incertezza e influenza nel modello, quindi è 
importante definire i loro valori nel processo di aggiornamento. Ci sono molti altri importanti 
parametri che possono essere considerati, come ad esempio la massa del materiale, ma la loro 
incertezza è relativamente bassa poiché la loro identificazione è stata effettuata in laboratorio. Ci sono 
diversi metodi per poter aggiornare il modello e in questo caso viene utilizzato quello proposto da 


Douglas e Reid. 


Calibration parameter Frequency optimized 
value 

Young's modulus of stone masonry Esm (MPa) 1560 

Young's modulus of brick masonry Epm (MPa) 1700 

Young's modulus of the reinforced masonry Em 1870 

(MPa) 

Lateral stiffness factor from Chapel 3.265 

Lateral stiffness factor from Sacristy 3.526 

Lateral stiffness factor from Tabernacle 3.063 


Tabella 37: Parametri per la calibrazione. 


La seconda fase di ottimizzazione si basa su un grande numero di modelli (tabella 38), nei quali il 
modulo di Young dei materiali varia entro un intervallo definito. I fattori di rigidità laterali non 


cambiano e le frequenze ottimizzate sono fisse. 
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Calibration parameter Min. Central Max. 


value value value 

Young's modulus of stone masonry Esm 780 1560 2330 
(MPa) 

Young's modulus of brick masonry Ebm 850 1700 2550 
(MPa) 

Young's modulus of the reinforced 940 1870 2810 


masonry Em (MPa) 


Tabella 38: Range (iniziale) di valori per i parametri, ottenuto con le frequenze del modello 


Il processo di aggiornamento del modello è stato fatto dal confronto dei parametri modali ottenuti 


sperimentalmente con quelli calcolati dalla perturbazione del modello agli elementi finiti. La distanza 


tra i modelli è stata calcolata con un’espressione, in cui oltre alle rispettive frequenze, compaiono dei 


fattori di ponderazione w. (figura 83). In questa seconda fase troviamo un range di valori ancora più 


ristretto (tabella 39). 


(a) Distance d considering w, = 0.50 yw, = 0.50 (b) Distance d considering w, = 0.50 y w,, = 0.50 
0.3 0.22 
NOMINAL FREQUENCY 
: 0.2 | vaLUES | OPTIMIZED VALUES 
0.25 NOMINAL FREQUENCY E, = 1390 MPa E, = 1560 MPa 
VALUES OPTIMIZED VALUES he ® 1690 MPa E * 1700 MPa 
E, 7 1390 MPa En ® 1560 MPa E, ® 1870 MPa E. = 1870 MPa 
Eon® 1690 MPa | Es, 1700 MPa 0.18 
0.2 [Em 1870 MPa E. = 1870 MPa 
0.16 
0.15 
0.14 
0.1 0.12 
0 100 200 300 400 500 600 700 800 900 300 320 340 360 380 400 420 440 460 480 500 
Model's number Model's number 


Figura 83: Distanza tra il modello sperimentale e numerico. (a) Distanza per tutti i modelli; (b) Zoom in un’area di interesse. 


Percent participation of frequency term in distance d 


100% 
90% 
Calibration parameter Min. Mac 80% 
value value 
Young'®s modulus of stone masonry EL (MPa) 1510 1550 70% 
Young®s modulus of brick masonry E, (MPa) 1540 1850 60% wi 
Young'®s modulus of the reinforced masonry Em 1780 2140 
(MPa) 50% —0.95 
— 0.90 
40% — 0.70 
è ; ; P 5 — 0,50 
Tabella 39: Range (finale) di valori per la calibrazione. » 
20% 
10% 
0% 
o 200 400 500 200 1000 
Madel's number 
Risultati e conclusioni Figura 84: Partecipazione del termine frequenza in base a w. 


Wo 

0.05. 
0.10 
0.30 
0.50 


Le informazioni ottenute dalle prove dinamiche in situ sono i dati di input per valutare la vulnerabilità 


sismica della struttura in esame, insieme ad un accurato modello agli elementi finiti. 
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I due step per l’aggiornamento del modello hanno prodotto risultati simili, anche se è importante 
notare che il costo del secondo è molto più grande. Il primo ottiene un unico valore per ogni parametro 
di calibrazione, mentre il secondo permette di ottenere intervalli di valori per ogni parametro. Questo 
può essere più coerente con la realtà della struttura, perché è più logico avere un intervallo valido per 
ciascun parametro di calibrazione, data la variabilità delle proprietà della struttura. Infine, sebbene 
non ci sia grande differenza tra i risultati dei due passaggi, vi è una maggiore certezza nell’ottenere i 
risultati nella seconda fase perché l’interazione raggiunta dall’utente con la regolazione del modello 


è molto maggiore. 


2.3.2.5 Chiesa di San Filippo Neri 


STATO Italia 
CITTÀ Macerata (MC) 
ANNO DI 
COSTRUZION 1697 d.C. 
E ; 
— i 
https://www.researchgate.net/publication/33549 | * ua 

FONTE È o È Figura 85: Vista dall'alto, Google Maps 


Tabella 40: Identificazione chiesa. 


Figura 86: Chiesa di San Filippo Neri 


Cenni storici e descrizione della struttura 


La chiesa è formata da un nartece rettangolare seguito dall’aula ovale principale, in cui si aprono 
quattro cappelle radiali, e una cupola ellittica con lanterna che la sormonta l’aula. 
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Nel 1697, durante la costruzione, le fondamenta 


originarie furono ampliate a causa della NI > 
- 
consistenza del suolo. Il 29 gennaio 1899 fu i 7 ia 


effettuato un sopralluogo nel seminterrato della 

casa adiacente, dove è stata scoperta una grotta con 

presenza di arenaria friabile inconsistente che si — 

estendeva sotto la chiesa. Anni più tardi venne 

costruita una volta in muratura all’interno della 

grotta e poi venne murata. Nel corso del 900 la 

chiesa subì diversi interventi, dovuto alla presenza Figura 87: Planimetria, stato fessurativo. 

di lesioni nella cupola, l’ultimo dei quali venne 

eseguito nel 2012, che comprendeva il rinforzo della cupola con materiali compositi, la riparazione 
di fessure e la posa in opera in un anello di rinforzo in acciaio nell’intradosso del tamburo. 

Dopo il terremoto Lazio-Marche-Umbria del 2016, la cupola della chiesa fu nuovamente danneggiata. 
Sono state condotte delle prove in situ, per la caratterizzazione meccanica dei materiali e per ricavare 


informazioni sul comportamento dinamico: 


e indagini endoscopiche; 

e tomografia sonica; 

e provacon martinetto piatto doppio; 

e prova con penetrometro dinamico sui giunti di malta; 


e provadi vibrazione ambientale. 


Vengono riportati solamente la prova con martinetto doppio e la prova di vibrazione ambientale, 


fondamentali per la calibrazione del modello FEM. 


Prova con martinetto piatto doppio 


Queste prove sono state eseguite al piano terra in modo 


da stimare il modulo elastico della muratura. In 


particolare, sono stati individuati tre diversi moduli Ue 
se e \ / gra # er LÉ, i 
elastici. E pl IF 


nz { 


Figura 88: Posizionamento dei martinetti doppi. 
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SUPERIOR CUT 


INFERIOR CUT 


4160,7 MPa 
7440,6 MPa 
M3 | 2896,1 MPa 


Figura 89: Esecuzione prova martinetti doppi. 


Prova di vibrazione ambientale 


Figura 90: Posizionamento degli accelerometri. 


L'analisi modale operativa è basata sull’assunzione delle vibrazioni 
casuali del rumore bianco. In questo scenario è sufficiente analizzare la 
risposta della struttura a vibrazioni sconosciute casuali (output-only 
modal identification) per estrapolare i parametri modali dell’edificio 
(frequenze naturali, vettori di forma modale, smorzamento modale). Il 
metodo usato per ricavare questi parametri è il metodo complesso nel 


dominio della frequenza dei minimi quadrati di poli riferimento. 


Per la visualizzazione dei risultati, viene considerato un modello 


geometrico (figura 91). Figura 91: Modello 
geometrico. 
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|_2,06$ Hz 
SF 


E=1,11% | Mode 2 E=2,84% | Mode 3 
_ Gr 


Global modes 


Local modes 


Figura 92: Modi sperimentali. 


Modellazione FEM 


Il modello FEM è basato su un accurato rilievo geometrico e rappresenta fedelmente tutti i volumi 
della chiesa. I solai in cemento armato dell’edificio annesso e il tetto orientale della chiesa sono 
modellati con i relativi vincoli, mentre l’influenza del complesso edilizio annesso è tenuto in 
considerazione con restraints orizzontali. Il modello considera un vincolo tridimensionale 
traslazionale per i nodi alla base, mentre i nodi laterali all’interfaccia edificio-terreno sono bloccati 
in direzione orizzontale. 

Il materiale è considerato omogeneo e isotropo con un comportamento elastico. Per differenziare i 


vari materiali 11 modello è stato diviso nelle principali macro-zone. Tabella 


Software Midas FEA 


Tipologia di Elemento |*Muratura: Elemento solido tridimensionale 
*Rinforzi in acciaio: elemento unidimensionale 


Dimensione e546049 elementi 3D 
€129 elementi 1D 
Calibratura Prova di vibrazione ambientale 


Tabella 41: Descrizione modello FEM. 
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Material Mass density [KN/m*] E [MPa] _v[-] 


Southern masonry — El 18 4160,7 0.21 
Northern masonry — E2 18 7440.6 0.30 
Eastern masonry — E3 18 2896.1 0.20 
Steel reinforcements 76.98 210000 0.3 


Tabella 42: Caratterizzazione dei materiali. 


Figura 93: Modello agli elementi finiti della chiesa. 


Model Updating 


Sulla base dei risultati sperimentali, è stata effettuata un’analisi di sensibilità al fine di abbinare 
l’analisi modale numerica con l’analisi modale operativa. Considerando la complessità dell’edificio 
e le incertezze legate alla tecnologia di costruzione, l’analisi si è incentrata sui modi globali della 
chiesa (primi sei modi di vibrazione). La procedura di aggiornamento del modello è basata sulla 


valutazione di due quantità: l'errore di frequenza relativa Dfî,j e il MAC (Modal Assurance Criterion). 


In questo studio si considera l’interazione con un terreno elastico (teoria di Winkler) e vengono 
considerate tre differenti gruppi di molle, con un range di ks tra 14400 e 36000 KN/m3. I risultati sono 
riportati in figura 94. 


Elastic Foundation Updated Model 
Subgrade reaction Variation MAC Modal Assurance Criterion Variation 
35 ] si 1,0 
to o... 2 0,8 
25 
î ed 0,6 
foro [8] mn 8 4 i 
error o 
15 0,4 
10 a °° > 
59 6 eris 02 i 
0 — rr]5 0,0 L 
150 200 250 1 2 3 4 5 6 7 8 
k. [Gpa] Modes 
dr mOricinal MUpdated 


Figura 94: A sinistra, risultati per la fondazione elastica in base a differenti ks ; a destra confronto tra modello 
originale e calibrato con MAC. 
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Frequency Variation 


Numerical frequencies [Hz] 


(0) 2 4 é 8 10 12 14 
Experimental frequencies [Hz] 


Figura 95: Confronto tra frequenze numeriche e 
sperimentali. 


Risultati e conclusioni 


I risultati dalle prove statiche hanno fornito importanti informazioni sui materiali e in particolare la 
prova con martinetti piatti doppi ha riportato una diversità del modulo elastico, cosa che è stata 
parzialmente validata dal modello FEM. La prova delle vibrazioni ambientali e la conseguente analisi 
modale operativa sono confermate come tecniche appropriate per la valutazione della risposta 
dinamica delle costruzioni storiche in muratura. 

La procedura di aggiornamento ha evidenziato la complessità del problema nella selezione ottimale 
dei parametri; l’applicazione di fondazioni elastiche è confermata la più rilevante per l’aggiornamento 


al fine di abbinare il comportamento modale numerico e sperimentale. 
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3. Capitolo 3 


Descrizione del caso studio: chiesa di Santa Maria in Via a Camerino 


3.1 Inquadramento storico 


La chiesa di Santa Maria in Via sorge nel centro storico di Camerino, in provincia di Macerata (Fig.96 
e 97). Le prime notizie della sua esistenza sono riportate nelle “Rationes Decimarum” del 1299 e in 
un testamento del 1341. Nel secolo XVI la chiesa era una costruzione povera ed irregolare, 
caratterizzata dalla presenza di una sacra icona della Madonna con Bambino, custodita in una piccola 
cappellina sorta come tante altre sulla strada dei pellegrinaggi romei. Successivamente nel XVII 
secolo il cardinale Angelo Giori acquistò le case limitrofe alla vecchia chiesa unendole alla canonica 
e oratorio preesistenti per rendere sufficientemente ampio il nuovo tempio, in modo da accogliere 
quanti venivano per venerare la sacra icona. 

L’attuale edificio fu costruito tra il 1639 e il 1642, dove la pianta irregolare ad andamento trapezoidale 
visibile all’esterno è dunque frutto dell’accorpamento di edifici preesistenti. Questo giustificherebbe 
la qualità delle murature il cui nucleo è particolarmente incoerente, mentre la fodera esterna venne 
aggiunta per regolarizzare 1 vari corpi di fabbrica. 

Fu consacrata nel 1654 dal Vescovo di Camerino monsignor Altieri alla presenza del Cardinale Giori 
che vi assistette in forma privata. Non si hanno notizie certe sulla costruzione del campanile della 
chiesa. Dall’osservazione delle caratteristiche architettoniche, probabilmente la sua realizzazione è 


contemporanea a quella della chiesa. 


Figura 96:L’antichissima città di Camerino nello Stato della Chiesa, 
Incisione di Thomas Salmon e pubblicata da Albrizzi a Venezia, 1757 


Figura 97: Archivio di Stato di Roma, Catasto Gregoriano, 1820. 
Chiesa parrocchiale di S. Maria in Via 
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Nel 1799 Camerino fu investito da un terremoto che provocò danni gravissimi alla città, come si può 
vedere in figura, e in particolar modo alle chiese. Una delle più danneggiate fu proprio Santa Maria 
in Via, dove crollò l’originaria volta in mattoni. Nei primi anni del XIX secolo la copertura crollata 


fu sostituita da un tetto a doppia falda sostenuto da capriate lignee e venne costruita anche una nuova 


cupola in camorcanna. 


Figura 98:Terremoto del 28 luglio 1799, distribuzione dei danni a Camerino. Da Capponi G., Il 
terremoto del 28 luglio 1799. Ricostruzione ed analisi dell'accaduto a Camerino. Tesi sperimentale di 
laurea in Sismologia storica, relatore Prof. Francesco Dramis, corr 


Nuovamente colpita dopo il sisma del 1997 fu oggetto di vari interventi di consolidamento statico e 


restauro architettonico da parte di varie imprese edili, per riaprire al culto nel settembre del 2006 
(Fig.99) 


Figura 99: Chiesa di Santa Maria in Via, prima del terremoto del 2016 
Nel 2016 a seguito del terremoto di agosto e ottobre venne chiuso tutto il centro storico di Camerino 


e ancora oggi rimane zona rossa. La chiesa in oggetto riportò gravi danni e attualmente è stata messa 


in sicurezza con un presidio in acciaio 
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3.2 Danni riportati a seguito del sisma 2016 


A seguito degli eventi sismici la chiesa subì diversi crolli e danni che la resero gravemente inagibile. 
Inoltre, nel gennaio del 2017 ci furono delle abbondanti nevicate che contribuirono al crollo di una 
buona parte del tiburio e della sovrastante copertura lasciando la chiesa esposta agli agenti atmosferici 
(Fig. 100); particolare preoccupazione era rivolta alle murature in pietra che si sono disgregate e che 


si sono trovate esposte al dilavamento di eventuali piogge battenti. 


Figura 100: Danni copertura. Vista esterna e interna. 


Il campanile è crollato alla seconda scossa del 26 ottobre alle ore 21.18, come testimonia l’orologio 
presente nella facciata che si è fermato qualche minuto dopo (Fig.101), e il basamento superstite si 
eleva per circa 4 metri sulla gronda del corpo principale ed è caratterizzato dalla presenza di 


incatenamenti e non è interessato da un quadro fessurativo. 


Figura 101: A sinistra campanile crollato, a destra particolare orologio. 
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La facciata presenta un meccanismo di ribaltamento (Fig.104) che interessa sia l’ordine superiore che 
inferiore e, dal quadro fessurativo e deformativo rilevato (Fig.102), si può collocare la cerniera appena 
sopra lo zoccolo del basamento, ad una quota circa +1,80 m dal calpestio del piazzale e una seconda 


cerniera in corrispondenza del cornicione del primo ordine. 


Figura 102: Risultato prova termografica 


La facciata nell’ordine inferiore ha subito una rotazione fuori dal piano, mentre il blocco superiore 
ha subito un moto rigido di traslazione. Questo meccanismo interessa anche una parte del tiburio, 


come testimoniano le lesioni presenti sullo stesso (Fig.103). 


Figura 103: Lesioni nel tiburio che testimoniano il meccanismo di ribaltamento 
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Figura 104: Meccanismo di ribaltamento della facciata 


Infatti, il distacco non si presenta solamente in corrispondenza della superficie di contatto tra il tiburio 
e il corpo di facciata, ma anche nel tiburio stesso dovuto alla presenza dei finestroni. Il paramento 
esterno in laterizio della facciata si è distaccato in diverse parti, come si può vedere in figura 105, 
dando origine a disgregazioni locali dovute al riempimento caotico della tipologia di muratura a sacco 
presente. Il fuori piombo è di circa 30 cm e un esame ravvicinato ha evidenziato una dislocazione del 
cornicione al primo ordine che sembra poggiare sul paramento in laterizio senza connessioni con la 


parete. 


Figura 105: Distacchi nel corpo di facciata. 
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Anche il portale in pietra arenaria ha subito dislocazioni molto significative con evidenti cedimenti 
degli elementi di architravatura. 

Il meccanismo di ribaltamento della facciata ha prodotto crolli molto significativi delle murature, 
soprattutto nei due lati del corpo di facciata, con formazione di lacune dovute sia al distacco del 
paramento in laterizio sia alla disgregazione della muratura di pietra. Inoltre, le pareti laterali sono 
interessate da un meccanismo secondario di ribaltamento fuori dal piano, come testimonia una lesione 
sulla parte destra della facciata. 

Oltre alle lesioni da distacco in corrispondenza della facciata e delle aperture più prossime alla stessa, 
i maggiori dissesti nel tiburio si hanno nella porzione posteriore che è crollata principalmente verso 
l’interno dell’aula per il danno progressivo innescatosi con gli eventi sismici e per le forti nevicate 
del gennaio 2017. Il crollo è da imputare sia alla pessima fattura della muratura e alla presenza del 
padiglione della copertura che presentava tre puntoni diagonali, sia al meccanismo di ribaltamento 
della facciata che ha coinvolto il tiburio, dove sono presenti dei tiranti che hanno quindi trascinato la 
parte posteriore. 

La rovina della parte di fondo del tiburio ha provocato il crollo parziale della porzione posteriore di 
copertura a padiglione lasciando molti elementi delle orditure e del manto in laterizio pericolanti. 


Oltre il 75 % della cupola in camorcanna è andata persa. 


Figura 106: Danni in copertura 
Le pareti interne della chiesa non sono caratterizzate da quadri fessurativi significativi, ad eccezione 
dei due altari laterali adiacenti al corpo di facciata che hanno subito danni importanti agli archi e alle 
volte. Il cornicione interno è stato danneggiato dalla caduta di blocchi dall’alto e presenta blocchi e 
porzioni di stucchi pericolanti. 
Le pareti laterali del corpo del corpo dell’aula e della parte absidale non sono interessate da crolli 0 


da quadri fessurativi particolarmente significativi, ad eccezione dello spigolo sud-ovest, al di sotto 
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del fusto del campanile, dove sono presenti delle lesioni diagonali, concentrate nelle fasce di piano 


tra un’apertura e l’altra. 


/ la E RZ: 
= o gira malk 
\ È O 
= i: mu] î 


sezione longitudinale AA' 


Figura 107: Prospetti e sezione, quadro fessurativo. 
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3.3. Intervento di messa in sicurezza 


L’obiettivo primario è quello di mettere in sicurezza le vie d’accesso e le aree pubbliche adiacenti 
alla chiesa con un intervento capace di ridurre al minimo l’occupazione degli spazi destinati alla 
viabilità. L’opera provvisionale è stata pensata per contrapporsi all’ulteriore mobilitazione dei 
cinematismi della facciata e del tiburio attivati nella sequenza sismica iniziata il 24 agosto 2016 e ai 
crolli che ne possano conseguire. L'intervento è organizzato in modo da controllare anche collassi 
parziali non facilmente prevedibili. 

Il danneggiamento manifestato e i meccanismi in atto costituiscono la prima fonte di informazione 
sulle vulnerabilità proprie della costruzione, sulla quale si basa l’intervento di messa in sicurezza. 
Viste le condizioni in cui versa l’edificio e la necessità di determinare in modo accurato le masse 
strutturali mobilitate sono stati condotti rilievi laser scanner e fotogrammetrici con l’ausilio di un 
drone. 

Gli interventi proposti consistono principalmente in un’operazione di incatenamento e trattenimento 
della scatola muraria del corpo di facciata e della porzione di tiburio più prossima alla stessa (Fig. 
108). Le dimensioni geometriche e le caratteristiche morfologiche delle pareti non permettono 
l’adozione delle soluzioni convenzionali solitamente usate per gli edifici multipiano ed è quindi stata 


adottata una soluzione originale, adeguata alle caratteristiche specifiche del manufatto. 


Figura 108: Funzionamento sistema messa in sicurezza 
L'intervento si completa con l’inserimento di un sistema tubi-giunti all’interno dell’aula e nella parte 
sommitale esterna del tiburio, che consenta la stabilizzazione e successiva ricostruzione della 
porzione di muratura crollata del tiburio, nonché la realizzazione di una copertura provvisoria, che 


impedisca ulteriori infiltrazioni di acque meteoriche all’interno dell’aula (Fig.109). In questo modo 
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la struttura realizzata sarà utilizzabile durante i successivi interventi di restauro e ha reso possibile la 


posa degli accelerometri nel cornicione per l’identificazione dinamica preliminare. 


Figura 109: Particolari messa in sicurezza, interno 


Per la facciata è stato predisposto un sistema esterno a graticcio (Fig.110), composto da tre ordini di 
elementi strutturali; il graticcio si collega alle estremità a due controventi concentrici in acciaio 
disposti lateralmente (prospetti Nord e Sud). 

Il sistema dei tralicci e dei controventi, fondato superficialmente su travi in calcestruzzo armato, è 


trattenuto da due ordini di tiranti longitudinali, rispettivamente n. 6 funi a trefolo @ 26 e n. 8 funi a 


Figura 110: Presidio in acciaio 
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trefolo 926, per ciascuno dei due lati, e due ordini di tiranti trasversali (n. 2 funi a trefolo @26 
ciascuno). I tiranti longitudinali corrono esternamente alle strutture della chiesa, percorrendone 
l’intero perimetro e trovando contrasto sul prospetto posteriore (prospetto Ovest) grazie ad un sistema 
di deviatori. I tiranti trasversali, invece, passano all’interno delle strutture murarie del tiburio e sono 
deputati a collegare mutuamente le ultime due elevazioni dei controventi reticolari, al fine di 
contrastare il possibile ribaltamento fuori dal piano delle porzioni laterali della facciata. 

Gli elementi principali che costituiscono il graticcio sono 4 travi reticolari tridimensionali in acciaio 
progettate affinché abbiano la necessaria resistenza e rigidezza con pesi propri contenuti. Esse si 
trovano sul fronte principale della chiesa e sono collegate all’estremità ai due controventi reticolari 
laterali. Delle quattro reticolari soltanto le prime due, partendo dal basso, sono collegate ai tiranti 
longitudinali. Gli elementi secondari sono costituiti da due livelli di profili verticali in acciaio, 
disposti sia in facciata, che lateralmente al di sopra e al di sotto del cornicione. Il terzo ordine degli 
elementi che compongono il graticcio è rappresentato da pannelli in legno tipo XLAM posti tra i 
profili verticali e la facciata (principale e laterale) della chiesa. Essi consentono il contrasto vero e 
proprio con le strutture murarie e sono adottati in quanto più facilmente adattabili alle anomalie 
geometriche delle superfici della muratura danneggiata. Il contrasto così operato vuole impedire la 
caduta di materiale in caso di nuove scosse ed ulteriori deformazioni prodotte dal danneggiamento 
sismico. 

Le azioni orizzontali che si sviluppano sulle pareti laterali e sulla parete posteriore, in conseguenza 
all’azione di tiro dei cavi longitudinali, potrebbero trovare dei punti di debolezza nelle nicchie e nelle 
aperture delle pareti che verranno rafforzate introducendo elementi in grado di assicurare la 


trasmissione delle forze orizzontali e presidi, quali sbadacchiature in legno. 


Figura 111: Messa in sicurezza facciata. 


76 


ELEMENT De PRESIDIO 
TO SINISTRO 
ino 


Î 
| 


ELEMENTI DI PRESIDIO 
POSTER 


Figura 112: Messa in sicurezza, prospetti. 
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3.4 Identificazione dinamica preliminare 


L’identificazione delle dinamiche globali della chiesa è stata effettuata attraverso prove di vibrazione 
ambientali impiegando 16 accelerometri e due configurazioni di sensori mostrate in figura 113a e 
113b. Per avere informazioni sul comportamento complessivo della chiesa, si è deciso di monitorare 
8 punti situati intorno all’ellisse del tiburio, a due altezze diverse dal livello del suolo. In ogni punto 
(da A a H) sono stati posizionati due accelerometri uniassiali con direzioni di misura lungo gli assi 
principali dell’aula della chiesa (chiamati direzioni x e y). In dettaglio, nella prima configurazione gli 
accelerometri sono posizionati in $ punti diversi (da 1A a 1H) a 12 m dal suolo, mentre nella seconda 
configurazione due accelerometri sono lasciati nella stessa posizione (1B), in modo da poter essere 
utilizzati come riferimento, mentre gli altri sono stati spostati a 19,7 m da terra; le posizioni 2C e 2H 


sono state escluse poiché inaccessibili. 


l n 2B 2D_2E 
2 2G 2F 
I +12,0 m | DA c=ass 
RES LE v — i i, 
i 1D, 1E FBI ir) IST 


Lr i; RIGA 
(b) 


Figura ai accelerometri: (a) Conf. 1: 8 stazioni di misura a +12m dal livello del terreno; (b) Conf. 2: stazione di 
riferimento e 6 stazioni a 19,7 m dal livello del terreno. 


Per quanto riguarda la strumentazione, sono stati utilizzati accelerometri PCB 393B31 con una 
sensibilità di 10V/g, moduli di conversione analogico-digitale NI 9234, un 9045 cRIO e tre 9185 
cDAQ. È stata creata una rete di sensori distribuiti posizionando un cRIO per la stazione di 
misurazione vicino al punto di riferimento 1B e tre cDAQ mobili vicino ai restanti punti di 
misurazione. La sincronizzazione è stata ottenuta tramite la tecnologia TSN che fornisce 
sincronizzazione temporale distribuita e comunicazione deterministica utilizzando reti Ethernet 
standard. 

Utilizzando i moduli NI 9234, i segnali analogici sono stati campionati a 2048 Hz. Prima 


dell’identificazione dei parametri dinamici, i segnali sono stati pre-elaborati rimuovendo la linea di 
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base, sottraendo il contributo ottenuto “adattando” il segnale con un polinomio di 3° grado, filtrando 
con un filtro passa-basso con frequenza di taglio di 25 Hz e ricampionando a una frequenza di 51,2 
Hz per ridurre la quantità di dati e rendere più rapida l’analisi successiva. 

L’identificazione dei parametri modali è stata effettuata utilizzando la tecnica Principal Component- 
Subspace Stochastic Identification (SSI-PC). Siccome le acquisizioni relative alle due configurazioni 
sono asincrone, per ottenere un corretto andamento del modo è stato necessario procedere con 
un’operazione di data merging; in questo caso sono riportati i risultati ottenuti con la tecnica POSER 
(Post Separate Estimation Re-scaling). 

La figura 114a il diagramma di stabilizzazione ottenuto dalle prove in Configurazione 1 e 
Configurazione 2. I modi stabili in termini di frequenza (differenza inferiore all’ 1%) e le forme modali 
(MAC maggiore del 95%) sono identificate da un cerchio nero pieno. Nella stessa figura sono riportati 
i diagrammi di frequenza-smorzamento, dove si possono apprezzare i valori dei rapporti di 


smorzamento relativi alle soluzioni stabili. 
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Figura 114: Risultati dell'’identificazione dinamica: (a) Diagrammi di stabilizzazione e trama Frequenza-Indice di Smorzamento; (b) 
parametri modali identificati. 
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La figura 114b mostra i parametri modali selezionati, cioè frequenze, rapporti di smorzamento e 
forme modali. A livello globale, sono state identificate 6 modi di vibrare. È interessante notare che, 
nonostante il complesso andamento delle fessure, la dinamica della struttura presenta i primi due modi 
traslazionali e ben disaccoppiati e un terzo modo torsionale, che ci si potrebbe aspettare per il sistema 


non danneggiato. Questo è chiaro anche dai valori della matrice autoMAC mostrati nella figura 115. 


2.13 4.8 0.0 4.3 
2.62 114 | 0.7 8.6 
3.33 |2% 


3.66 | 4.8 11.4 


3.96 | 0.0 0.7 2.0 


4.36 | 4.3 8.6 0.5 1.1 


f[Hz] 2.13. 2.62 3.33 3.66 3.96 4.36 
Figura 115: Matrice autoMAC 


Per quanto riguarda i modi superiori, questi necessitano di ulteriori approfondimenti ma nel 


complesso sembrano coinvolgere le dinamiche locali del tiburio e della facciata. 
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3.5 Sistema di monitoraggio 


La chiesa di Santa Maria in Via è stata severamente danneggiata a seguito del terremoto del 2016 che 
ha colpito il centro Italia, quindi è stata messa in sicurezza ed è stato previsto un monitoraggio per 
controllare l'evoluzione del danno. 
Attualmente è stata predisposta solamente una parte del sistema di monitoraggio per studiare la 
dinamica della facciata (Fig. 116), ovvero: 

- 2 accelerometri piezoelettrici PCB 393A03 con una sensibilità di 1V/g e la direzione di misura 

perpendicolare alla facciata a +12.6 m di altezza dal terreno; 
- una termocoppia di tipo T posizionata su un elemento del sistema di sicurezza in acciaio e una 


termocoppia dentro la chiesa. 


© PCB 393408 ® Thermocouple 


Figura 116: Layout dei primi sensori de sistema di monitoraggio 


Il software di acquisizione registra continuamente i dati, in tempo reale, e li invia tramite il protocollo 
FTP al server collocato alla facoltà di Ingegneria in Ancona. Questo approccio, grazie allo streaming 
continuo dei dati e al determinismo garantito dal sistema di acquisizione, permette di indagare in 
breve tempo eventuali dipendenze dalle condizioni ambientali, ma anche di registrare la risposta ad 
eventi sismici senza il rischio di perdere campioni di dati. 

I dati vengono organizzati in file da 20 minuti e i parametri modali vengono identificati e selezionati 
automaticamente. L’algoritmo di identificazione automatica prevede una prima fase in cui si ottiene 
il diagramma di stabilizzazione applicando l’algoritmo SSI-PC, dopodiché i parametri modali 
vengono selezionati mediante un algoritmo basato su un AHC (Agglomerative Hierarchical 
Clustering). La distanza tra le soluzioni del diagramma di stabilizzazione è calcolata utilizzando 


l’equazione proposta da Magalhaes. 
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Una chiara dipendenza delle frequenze dalle temperature non è ancora chiara, probabilmente a causa 


dei pochi dati a disposizione, poiché i dati elaborati di questa prima parte del monitoraggio fanno 


riferimento solamente a un mese. Inoltre, siccome la chiesa si trova in zona rossa, non è soggetta a 


molte vibrazioni ambientali e quindi l’identificazione automatica non riesce sempre a identificare 


tutti i modi che contribuiscono all’identificazione dinamica di tutta la struttura. 


Nei prossimi mesi verrà completato il sistema di monitoraggio con i sensori rimanenti, che sono: 


qualche accelerometro nel tiburio; 

qualche accelerometro sulla facciata; 

2 velocimetri a livello di fondazione (fondamentali in caso di eventi sismici); 

accelerometri wireless mems per monitorare i cavi a trefoli di acciaio utilizzati nel sistema di 
messa in sicurezza; 

sensori di temperatura e umidità all’interno e all’esterno, poiché le condizioni metereologiche 
possono influenzare la dinamica della struttura; 

inclinometri e potenziometri lineari, per monitorare l’evoluzione della cinematica della 


facciata. 
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4 Capitolo 4 
Modellazione FEM e Model Updating 


4.1 Rilievo geometrico e descrizione della struttura 


Per rappresentare al meglio la geometria, gli spazi e quindi le masse effettivamente presenti 
all’interno della struttura, è stato eseguito un rilievo geometrico, grazie al quale è stato possibile 
costruire un modello tridimensionale dettagliato (Fig. 117-118). 
La chiesa è costituita da un unico corpo di fabbrica dal quale si elevano il tiburio, a pianta ottagonale 
irregolare, 11 campanile, crollato a seguito del terremoto del 26 ottobre 2016, e il corpo della facciata. 
L’aula presenta una pianta centrale a forma ellittica, in cui l’asse maggiore di lunghezza 28 m 
comprende l’ingresso e l’altare maggiore caratterizzato da un coro profondo absidato. Lateralmente 
si innestano quattro cappelle radiali a nicchia emiciclica in spessore, che lasciano il pieno murario 
sull’asse minore, di lunghezza 15 m, dove sono presenti due coretti, al di sopra dei due accessi laterali. 
L’ingresso è compreso da due locali, il battistero a sinistra e la Cappella del Crocifisso a destra, 
mentre ai lati dell’abside sono presenti l’oratorio della confraternita, purtroppo inaccessibile, e la 
sacrestia, dalla pianta ottagonale, dalla quale si accede all’edicola della Vergine, ovvero una delle 
quattro cappelle radiali, aperta sulla chiesa. 
All’esterno la Chiesa, nelle sue forme architettoniche estremamente semplici e lineari, non lascia 


intuire la complessa composizione dello spazio interno. Essendo frutto di accorpamenti di edifici nel 
corso degli anni, la muratura del corpo della chiesa è prevalentemente a sacco, con uno spessore non 


uniforme in tutto l’edificio, e presenta una pianta trapezoidale. 
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Figura 117: Planimetria a -1.00m a sinistra, a +1.50 ma destra, con riferimento foto 
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Figura 118: Planimetria a +8.65 m, a + 12.00m e a 19.70 m, con riferimento foto 


Il corpo di facciata è realizzato in murature differenti per qualità e dimensioni dei mattoni e presenta 


una pianta trapezoidale isoscele con lunghezza media di circa 16 m ed altezza di circa 6 m. È 


organizzata su due piani sovrapposti (due ordini) e divisa da quattro paraste in tre campate, con la 


centrale più risaltata e con un frontone triangolare in cima (Fig.119). L’ordine inferiore è inglobato 


nella pianta trapezoidale della chiesa mentre l’ordine superiore si eleva, in adiacenza al tiburio, 


rispetto alla copertura degli altari laterali per un’altezza massima di circa 23 m. Al suo interno il corpo 


di facciata presenta tre orizzontamenti; il primo, a circa +6.60 m dal calpestio dell’aula, è costituito 


da volte a botte con generatrice parallela alla facciata, il secondo, a circa +12 m dal calpestio dell’aula 


e il terzo alla quota dell’orologio, ovvero circa +17 m dal calpestio dell’aula. In pianta è suddiviso in 


tre vani, a ciascuna delle elevazioni, separati da due pareti di spina in muratura caratterizzate da porte 


allineate sulla verticale. All’interno della parete di facciata, sul lato sinistro guardando dall’esterno, 
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risulta inglobata una scala a chiocciola che permette di accedere alla cantoria e al piano superiore. 


Questa scala produce una discontinuità significativa all’interno della muratura della facciata. 


+17.00m 


Figura 119: Quota orizzontamenti corpo di facciata 


Il tiburio si eleva per oltre 8 metri dal corpo principale della chiesa ed è caratterizzato dalla presenza 
di quattro finestroni in corrispondenza degli altari laterali. Esternamente si possono notare dei 
contrafforti localizzati negli angoli per un totale di 6. 

Il campanile crollato a seguito del sisma presentava una pianta ottagonale e cella campanaria piuttosto 
slanciata. 

Dalla sacrestia si accede ad una scala a chiocciola in pietra, situata sotto il campanile, dalla quale si 
arriva al piano superiore a circa +8.65 m dal piano di calpestio, dove è presente l’appartamento privato 
del sacerdote. La scala arriva fino a -2.77 m, dove è situata un’entrata posteriore. 

L’appartamento si trova nella parte posteriore della struttura e si sviluppa in due piani. Al primo 
piano, proprio sopra la sacrestia, sono collocati la cucina e il soggiorno, dal quale si arriva alla zona 
notte tramite un passaggio dietro l’abside dell’altare maggiore. Dal soggiorno si arriva anche al piano 
superiore dove sono presenti altre camere da letto e servizi igienici, che si trovano sopra la zona notte 
del primo piano e sopra l’altare maggiore. L'appartamento è collegato al corpo di facciata tramite un 
passaggio che dalla cucina arriva alla seconda elevazione della facciata, ovvero a +12.00 m. Allo 
stesso modo, dalla parte opposta della chiesa, si arriva alla stessa quota della facciata passando dal 
secondo piano dell’appartamento, con la differenza che non vi è un’apertura che giunge direttamente 
alla facciata, ovvero il passaggio è senza uscita. La scala a chiocciola non è l’unico collegamento 


verticale; dalla sacrestia si può giungere all’appartamento anche tramite un ascensore. 
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Particolare scala a chiocciola 
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PROSPETTO ANTERIORE PROSPETTO POSTERIORE 
NORD-EST SUD-OVEST 


PROSPETTO LATERALE PROSPETTO LATERALE 
NORD-OVEST SUD-EST 


Figura 120: Prospetti 


Grazie al rilievo geometrico è stato possibile avere tutte le informazioni e le corrette misure per 
disegnare il modello tridimensionale, anche se la struttura presenta una geometria molto complessa e 
particolare. Inoltre, alcune parti inaccessibili sono state ipotizzate facendo riferimento al progetto 
dell’ultima ristrutturazione dove purtroppo alcune informazioni sono mancanti. 

Dopo aver riportato le piante alle diverse quote, sono stati creati i volumi delle varie parti della 


struttura, in modo da sottrarli al 3D pieno ed ottenere i vuoti (Fig.121-122). 


Figura 121:Volumi pieni e struttura con i volumi tolti. 
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Figura 122: Volumi pieni e struttura con i volumi tolti. 


Figura 123: Modello trasparente, vista prospetto e pianta 
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4.2 Rilievo materico 


4.2.1 Murature 


La Chiesa presenta un rivestimento esterno continuo in laterizio, realizzato in modo da rendere 
omogenea la struttura visti i continui accorpamenti negli anni. Non sono state effettuate prove 
meccaniche in situ o in laboratorio, ma un’accurata indagine visiva ha permesso l’individuazione 
delle diverse murature. 

Sono presenti quattro tipologie di muratura diverse sia per materiali (pietra arenaria e pietra dura) che 
per forma (Fig.123). Per tutte e quattro le tipologie di muratura è stato assunto un peso specifico di 


19 KN/mì, valore indicato nel progetto della più recente ristrutturazione. 


Muratura 1 


Muratura 2 


Muratura 3 


Muratura 4 


Figura 124: Tipologie di muratura presenti 


Stima modulo elastico 


Il modulo elastico delle quattro tipologie di muratura è stato stimato in base ai valori presenti nella 
tabella C8.5.I nella circolare applicativa alle NTC18 (Fig.125), data la mancanza delle prove in situ. 
Per le murature 2 e 3 sono stati presi direttamente i valori medi, visto che appartenevano ad una delle 
categorie in tabella, invece per le restanti murature si è proceduto diversamente poiché sono murature 
composte da due paramenti, con l’eventuale presenza di nucleo interno. Sono stati presi i valori medi 
dei due paramenti ed è stata calcolata nuovamente la media tra i due. In caso di nucleo interno, il 


modulo elastico viene moltiplicato per 0,8. 
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f 
Tipologia di muratura (N/mm?) 


Muratura in pietrame disordinata (dottoli, pietre erratiche e 10-20 0,018-0,032 IC 690-1050 230-350 
irregolari) 

Muratura a conci sbozzati, con paramenti di spessore 0,035-0,051 1020-1440 340 
disomogeneo (*) 

Muratura in pietre a spacco con buona tessitura 2,6-38 0,056-0,074 a 1500-1980 


Muratura irregolare di pietra tenera (tufo, calcarenite, ecc.,) 1,4-2,2 0,028-0,042 EE 300-420 
13+16(**) 
.» a conci regolari di pietra tenera (tufo, calcarenite, 2,032 0,04-0,08 [010019 | 1200-1620 


Muratura a blocchi lapidei squadrati 5,882 0,09-0,12 | 0,18-0,28 2 22 
Muratura in mattoni pieni e malta di calce (***) 2,643 0,05-0,13 |0,13027 | 1200-1800 | 400-600 | 


Muratura in mattoni semipieni con malta cementizia 
5,0-8,0 0,08-0,17 0,20-0,36 875-1400 
(es: doppio UNI foratura <40%) 


Figura 125: Tabella C8.5.I, Circolare applicativa alle norme tecniche per le costruzioni 
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Muratura 1 


Muratura a sacco, composta da un paramento esterno di mattoni in laterizio di spessore 1 o 2 teste e 
un paramento interno di muratura a conci sbozzati, con riempimento caotico. Questa tipologia di 
muratura è presente nella parte esterna del corpo di facciata. Dalle immagini si può notare che non è 
una muratura di buona qualità, soprattutto per la presenza di un riempimento caotico, nel quale è 
possibile individuare una grande quantità di malta e materiale di scarto. In figura si nota anche la 
presenza di coppi in laterizio usati in copertura. Il paramento esterno presenta dei distaccamenti, 


dovuti all’assente connessione con il paramento interno. 


Figura 126: Muratura 1 
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Modulo elastico El 


paramenti di spessore 
disomogeneo 


2 
; E (N/mm?) E (N/mm?) DT, 
Paramenti . 4 Valore 
min-max Valore medio . 
medio*0,8 

Muratura In mattoni pieni in 1200-1800 1500 

laterizio e malta di calce 

Muratura a conci sbozzati con 1092 
1020-1440 1230 


Tabella 42: Modulo elastico El 
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Muratura 2 


Muratura a conci sbozzati con paramenti di spessore disomogeneo. Questa tipologia di muratura è 


presente nella parte interna del corpo di facciata e nei due muri di collegamento tra la parte interna e 


la parte esterna, e nelle pareti dei locali privati della chiesa. 


Figura 128: Muratura 2 


Modulo elastico E2 


paramenti di spessore 
disomogeneo 


2 2 
Paramenti E (N/mm?) E2 (N/mm°) 
min-max 
Muratura a conci sbozzati con 
1020-1440 1230 


Tabella 43: Modulo elastico E2 


Figura 129: Immagini muratura 2 


Muratura 3 


Muratura a sacco con paramenti in laterizio con il nucleo interno composto da muratura a conci 
sbozzati e riempimento caotico. La tipologia di muratura 4 è la muratura del corpo della chiesa ed ha 
uno spessore molto variabile, mantenendo costante il paramento esterno in laterizio visto che funge 


principalmente da rivestimento. 


Figura 130: Muratura 3 


Figura 131: Immagini muratura 3 
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Modulo elastico E3 


2 E (N/mm?) E3 (N/mm?) 
Paramenti E (N/mm ) Valore medio Valore 
min-max . 
medio*0,8 
I e 1200-1800 1500 1200 
laterizio e malta di calce 


Muratura 4 


Muratura con paramento interno in laterizio (spessore di circa 1 testa) e paramento esterno in muratura 
a conci sbozzati con presenza di listature. La muratura è presente nelle strutture verticali del tiburio. 
In alcune parti si notano delle discontinuità, frutto di varie modifiche nel corso degli anni. A 


differenza delle altre tipologie, in questo caso il paramento esterno presenta delle listature visibili in 


Tabella 44: Modulo elastico E3 


figura 133. 
Figura 132: Muratura 4 
Modulo elastico E4 
. E (N/mm?) E (N/mm?) E4 (N/mm?) 
PARI min-max Valore medio Valore medio 
Muratura In mattoni pieni in 1200-1800 1500 
laterizio e malta di calce 
Muratura a conci sbozzati con 1365 
paramenti di spessore 1020-1440 1230 
disomogeneo 


Tabella 45: Modulo elastico E4 
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xe RAPE BRE 


Figura 133: ii muratura 4 


MURATURA MODULO LL SINICO E 
(N/mm?) 
1 1092 
2 1230 
3 1200 
4 1365 


Tabella 46: Moduli elastici utilizzati per Modello Numerico I 
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4.2.2 Solai in legno 


All’interno della struttura sono presenti degli orizzontamenti composti da solai in legno. Prendendo 


come riferimento i progetti dell’ultimo intervento realizzato (Fig.134), è stato possibile individuare 


con precisione la presenza del solaio e le sue dimensioni e composizione. I solai in legno si trovano 


negli orizzontamenti della seconda e terza elevazione della facciata e in tutti gli orizzontamenti 


dell’appartamento privato presente nella parte posteriore della struttura. 


cena I \SOLAIO:S0 | |sOLAIOSI — 


pavimento 


roassetto alle: 
perso (em #) 
pl 


| filetti Coma Rx4) 
cavi 


Scala 1:20 


| È FIN — = ner: 
i | î SOLAIO S2 SOLAIO S3 
xi perte pavimento 


piane 


| 
masso all: | 

[ | DER —r —————m | 
pismlle | x 

1) | Eletti (ona 814) 


ALL LIT[[onas, x A î Tia | Alcii (ma 8x0) __ ge” fe 
3 cavi 


| Scala 1:20 Salle 120 


Figura 134: Tavole intervento di rifacimento solai in legno, 1997. 


Mmissetto alleg- _ 
gerito (om) _——— 


Il materiale è stato assunto con un comportamento omogeneo e isotropo, con un modulo elastico pari 


a 600 N/mm? e un peso specifico di 800 kg/m?. 


Figura 135: Distribuzione solai in legno nella chiesa 
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4.2.3 Copertura 


La copertura del tiburio è a padiglione ed è sostenuta da 3 capriate in legno e 4 capriate in acciaio. 
Prendendo come riferimento una tavola del progetto dell’ultimo intervento in copertura (Fig.137), è 
stata eseguita una stima dei pesi propri delle due tipologie di capriate in base alle loro dimensioni, 
ipotizzando un peso specifico delle travi in legno di 800 kg/m}. Il peso di progetto del manto di 
copertura è di 510 kg/m? (Fig.136). Si riporta il calcolo (metodo semplificato) delle forze puntuali 
all’appoggio delle capriate da applicare al modello 1. 


RELAZIONE TECNICA IN RIFERIMENTO ALLE DIRETTIVE 
TECNICHE DI CUI ALLA DELIBERAZIONE DELLA 
G.R. N° 78/99 PR/CBC DEL 18/1/1999 ALL.A/1. 
CHIESA DI SANTA MARIA IN VIA 


Forza sismica: F=FCxRxexBxIx W 


R=1 

Camerino: e=F,=1.2 
p=4 

I=1.2 


F=Cx1x1.2x4x1.2x W=Cx5.76 x W 
y muratura= 1600 Kg/m' 
ox=3 Kg/em? 
Carichi del tetto: stato attuale 
Neve = 240Kg/mq (Camerino = 657m s.l.m.; Qsk=160 + 300(657-200)/1000= 
=300Kg/mq: per un coefficiente dato dalla pendenza delle falde di circa 16° = 0.80 
otteniamo 300Kg/mq x 0.80= 240Kg/mq) 
Orditura lignea= 30Kg/mq 
Pianelle= 60Kg/mq 
ci oppi = 80Kg/mq 
È 40 


DINLTPRO guasta cimqure ni Leno pr na è co coo 


a 137: Tavola di progetto intervento capriate, 1997 
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Figura 139: Particolare appoggio delle due capriate 


. . . Figura 138: Identificazione in pianta capriate in legno e in acciaio 
Peso proprio capriate in legno 1, 2, 3: S Ron + 9 


Tabella 47: Calcolo peso proprio capriata I 


Puntone 0,3 0,3 13,16 | 800 | 947,52 

MIGHACO:E 0,3 0,3 6,48 800 | 466,56 

saettoni 

Tirante 0,36 0,36 | 13,37 | 800 |1386,20 
TOT. |2800,28 
TOT. /2 | 1400,14 


Tabella 48: Calcolo peso proprio capriata 2 


Puntone 0,3 0,3 14,66 | 800 |1055,52 

MGNACO:e 0,3 0,3 6,39 800 | 460,08 

saettoni 

Tirante 0,41 0,41 15 800 |2017,20 
TOT. |3532,80 
TOT. /2 | 1766,40 


Tabella 49: Calcolo peso proprio capriata 3 


Pimione 0,3 03 | 1298 | 800 | 934,56 

MONACO 0,3 03 | 642 | 800 | 462,24 

saettoni 

rTiranfe 0,36 | 0,36 | 13,29 | 800 |137791 
TOT. |2774,71 
TOT. /2 | 1387,35 


- Peso proprio capriate in acciaio 4, 5, 6, 7: 


Tabelle 50-51: Calcolo peso proprio capriata 4-6 


Tabelle 52-53: Calcolo peso proprio capriata 5-7 


HEA180 13,16 93 467,18 
HEA120 13,3% 19,9 266,06 
HEA100 6,48 16,7 108,22 
TOT. 841,46 
TOT./2 | 420,73 
HEA180 15,32 35,5 943,93 
HEA120 14,02 19,9 278,92 
HEA100 9,93 16,7 165,90 
TOT. 988,75 
TOT./2 | 494,37 


- Forze puntuali agli appoggi. 


HEA180 15,32 SS 543,93 
HEA120 14,02 19,9 278,92 
HEA100 9,93 16,7 165,90 
TOT. 988,75 
TOT./2 | 494,37 
HEA180 13,43 sole 467,62 
HEA120 12,22 19,9 243,10 
HEA100 8,966 16,7 149,73 
TOT. 869,45 
TOT./2 | 434,73 


Le forze puntuali vengono calcolate secondo un approccio semplificato: si considera un interasse di 


2,40 me si somma il contributo del peso proprio della capriata con il peso puntuale della copertura. 


Visto che la capriata 1 e 4 sono praticamente affiancate, viene considerato un unico punto di appoggio, 


sommando entrambi i pesi propri. 


1+4 > [(510*2,40) *13,37] /2 + 1400,14 + 420,73 = 10003 kg = 98,1 KN 
2 > [(510*2,40) #15] /2 + 1766,4 = 10946 kg = 107,34 KN 

3 > [(510*2,40) *13,29] /2 + 1387,35 = 9520,8 kg = 95,37 KN 
5 —> [(510*2,40) *14,02] /2 + 494,37 = 9074,6 kg = 88,99 KN 
6 > [(510*2,40) *14,02] /2 + 494,37 = 9074,6 kg = 88,99 KN 
7 [(510*2,40) *12,22] /2 + 434,73 = 7913,4 kg = 77,6 kN 


4.2.4 Regioni ripristinate 


Per proteggere l’interno e la muratura dagli agenti atmosferici, alcune regioni sono state ripristinate 
con un intervento di messa in sicurezza, che consiste in una struttura reticolare di tubolari in acciaio 
(Fig.140). 

Siccome potrebbero influire nel comportamento dinamico della chiesa, queste regioni sono state 
considerate nel modello, ipotizzando un comportamento più rigido rispetto alla muratura. 


Il valore del modulo è stato assunto pari a 1500 N/mm?, con un peso specifico di 1000 kg/m}. 


Figura 140: Identificazione regioni ripristinate con relative foto 
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4.3 Descrizione modello FEM 


Una volta completato il modello 3D in Autocad, è stato possibile sviluppare il modello numerico nel 
software Ansys. Per prima cosa è stato importato in Ansys Workbench, per essere corretto e per 
eliminare tutte le piccole imperfezioni per consentire l’esecuzione di una mesh precisa. 

Nel modello di partenza (Fig.141) viene considerata la copertura in termini di massa e non in termini 
di rigidezza e si ipotizza che la soletta dei solai in legno sia collaborante, visto che sono stati rifatti e 
sistemati nell’ultimo intervento a seguito del terremoto del 97. A seguito di questa ipotesi i solai in 
legno sono stati modellati come elementi shell, con uno spessore equivalente di 30 cm. Queste due 
assunzioni verranno poi studiate e valutate, con la costruzione di altri due modelli. 

Il modello è stato diviso in diversi blocchi a seconda dalla tipologia di muratura e inoltre, il blocco 
rappresentante il corpo della chiesa, viene suddiviso in tanti corpi vista la sua complessità. In questo 
modo la generazione della mesh è risultata più semplice in quanto il modello risulta costituito da 


blocchi di geometria più regolare. 


ANSYS 
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pia Y 
0.000 10.000 20.000 (m) x 


Figura 141: Modello agli elementi finiti della chiesa di Santa Maria in Via. Modello MN 


Il modello numerico è quindi composto da 535850 nodi e 308823 elementi con una mesh di 
dimensione 0.5 m, poiché una mesh più fine produce una risposta più precisa, anche se incrementa il 
tempo di calcolo e la memoria richiesta. In particolare, viene utilizzato un elemento tetraedrico a 4 


nodi. 
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I nodi alla base sono incastrati a terra e i materiali utilizzati sono considerati omogenei e isotropi con 
un comportamento elastico lineare. 
Per la caratterizzazione meccanica dei materiali non sono state eseguite delle prove e le varie tipologie 
della muratura sono state individuate con un’indagine visiva. Quindi, ci sono delle incertezze relative 
ai materiali e le caratteristiche elastiche sono state valutate prendendo i valori medi delle tipologie 
simili presenti in letteratura. Il modello risultante viene chiamato Modello Numerico (NM). 
Nel modello numerico sono stati inseriti 6 parametri diversi a seconda della tipologia di muratura o 
elemento costruttivo, e sono i valori su cui si concentrerà l’aggiornamento del modello. 
I parametri sono: 

- Modulo elastico muratura 1, E]; 

- Modulo elastico muratura 2, E2; 

- Modulo elastico muratura 3, E3; 

- Modulo elastico muratura 4, E4; 

- Modulo elastico regioni ripristinate, ES; 


- Modulo elastico solai in legno, ES. 


Dopo aver eseguito varie analisi, si è scelto di valutare anche l’effetto della rigidezza dato dalle 
capriate di copertura, inserendo nel modello 6 elementi rigidi senza massa, visto che i carichi di 
copertura sono già stati assegnati. Questi elementi sono posizionati allo stesso modo delle capriate, 
con il rilascio alle estremità, in modo da simulare il vincolo di appoggio delle travi. Il modello 


risultante viene chiamato Modello Numerico Copertura (NMC) (Fig.142). 
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0.000 10.000 20.000 (m) 
5,000 15.000 
Figura 142: Modello agli elementi finiti della chiesa di Santa Maria in Via. Modello MNC 
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Infine, si è scelto di andare a valutare l’ipotesi della soletta collaborante, costruendo un modello in 
cui vengono modellate le travi in legno del solaio, al posto dell’elemento shell. Il modello risultante 
viene chiamato Modello Numerico Trave (NMTY), nel quale vengono inseriti anche gli elementi rigidi 
fungenti da travi di copertura (Fig. 143). 

La muratura viene modellata con la stessa mesh di 0.5 m, mentre le travi vengono modellate con una 


mesh più fine, di 0.1 m. Il modello sarà composto da 887201 nodi e 358384 elementi. 
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7.500 22.500 
Figura 143: Modello agli elementi finiti della chiesa di Santa Maria in Via. Modello MNT 
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0.000 15.000 30.000 (mm) pix 
rr P I __—_—__ ni Y 


Figura 144: Modello MNT 
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4.4 Modello 1 


Il modello 1 (MNI1) prende in considerazione i moduli elastici valutati in fase d’indagine visiva e 


considerando il materiale con un comportamento omogeneo e isotropo. 


Figura 145: Identificazione materiali 


Materiale Peso specifico [KN/mÎ] E [N/mm?] Vv 
Muratura 1 19 1092 0.25 
Muratura 2 19 1230 0.25 
Muratura 3 19 1200 0.25 
Muratura 4 19 1365 0.25 
Parti ripristinate 10 1500 0.25 
Solai legno 8 600 0.25 


Tabella 54: Caratterizzazione dei materiali. Moduli elastici di partenza. 


Come prima cosa è stata svolta un’analisi statica in modo da poter inserire i carichi da copertura e 
poterli considerare poi nell’analisi modale (Fig.146). Infatti, il codice Ansys permette di creare dei 
collegamenti tra le due tipologie di analisi, visto che nell’analisi modale non si possono inserire i 
carichi esterni. 

Una volta eseguita l’analisi statica, si è proceduto con l’analisi modale per andare a determinare i 
modi di vibrare della struttura. L’analisi è stata eseguita con 15, 50, 80 e 100 modi per cercare di 
mobilitare più massa possibile. In tutti i casi la massa totale partecipante non cambiava di molto e 
anche con 100 modi non si è raggiunto nemmeno 1°80%. Questo è dovuto al tipo di struttura e al 


corpo della chiesa, in cui è presente molta massa difficile da coinvolgere e mobilitare (Fig.147). 
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Force: 98100N 
Force 6: 95370N 
Force 5: 95370N 

5 Force 7: 88990N 

{NN Force 9: 88990N 

[EN Force 10: 89990N 


Figura 146: Assegnazione del peso proprio e dei carichi da copertura, a sinistra, risultato analisi statica, a destra. 


EFFECTIVE MASS 


A x-DIR RATIO$ Y-DIR RATIO% Z-DIR RATIO% 
1° Modo 1238, 0.02 0.1663E+07 22.34 114.8 0.00 
2° Modo) 0.1011E+07 13.58 4193. 0.06 6537. 0.09 
3° Modo | 0.1362E+06 1.83 0.2415E+07 32.44 0.4244 0.00 

TOT 0.5729E+07 76.97 0.5806E+07 78.00 0.4155E+07 55.82 


Figura 147: Massa partecipante 


I risultati dell’analisi dinamica si possono vedere nelle figure 148-149, dove sono illustrati i primi tre 
modi di vibrare, ovvero i due traslazionali in X e Y e il torsionale in Z, e altri due modi, uno 
traslazionale in Y e l’altro torsionale in Z, ricavati poiché somiglianti al 5° e 6° modo del modello 
sperimentale. La tabella conferma la difficoltà di mobilitare la massa, infatti si può notare come i 
modi coinvolgano principalmente la facciata e il tiburio. Le direzioni che sono indicate nei modi di 
vibrare non corrispondono alle direzioni in tabella poiché, per essere coerenti con i risultati delle 


prove sperimentali, la X e la Y sono state scambiate. 


1° Modo 2° Modo 2° Modo 
Trasl. X Trasl. Y Tors. Z 
F3.05 Hz F3.33 Hz FA4.08 Hz 


Figura 148: Forme modali e frequenze dei primi tre modi di vibrare 
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O N 


8° Modo 9° Modo 
Trasl. Y Tors. Z 
5.84 Hz fF6.28 Hz 


Figura 149: Forme modali e frequenze dell'ottavo e nono modo di vibrare. 
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4.5 Confronto con 1 dati sperimentali 


Sono state eseguite delle prove alle vibrazioni ambientali (Fig.150) per poter effettuare una 
preliminare identificazione dinamica. In totale sono stati individuati 6 modi di vibrare e dai risultati 
del modello numerico, risulta una corrispondenza di 5 modi di vibrare del modello sperimentale. In 
particolare, la corrispondenza vi è nei primi tre modi di vibrare di entrambi i modelli, e tra il 5 e 1°8 


e tra il 6 e il 9 del modello sperimentale e numerico, rispettivamente. 


Configurazione 1 Configurazione 2 
+12.00m +17.00m 


Figura 150: Posizionamento accelerometri, punti di misura 


1° Modo 2° Modo 2° Modo 5° Modo 6° Modo 
Trasl. X Trasl. Y Tors. Z Trasl. Y Tors. Z 
f2.13 Hz F2.62 Hz {3.33 Hz 3.96 Hz F4.36 Hz 


Figura 151: Forme modali e frequenze dei cinque modi di vibrare del modello sperimentale 
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Modo fis [Hz] fwni[Hz] Afi [%] 


1 2.13 3.05 -30 
2 2.62 3,33 -21 
3 3.33 4.08 -18 
4 3.96 5.84 -37 
5 4.36 6.28 -25 


Tabella 55: Confronto delle frequenze tra il modello sperimentale (MS) e il modello numerico (MNI) 


Dalla tabella si può notare che il modello numerico (MN) è più rigido del modello sperimentale (MS), 
poiché le frequenze sono maggiori. Per valutare la distanza tra i due modelli è stato calcolato l’indice 
MAC (Modal Assurance Criterion), dove si può notare una discreta corrispondenza per i primi 3 
modi, mentre per gli altri due modi i valori sono molto inferiori. Per questo motivo, l’analisi si è 
concentrata solamente nei i primi 3. 

Nel modello numerico sono stati individuati i nodi corrispondenti alla posizione degli accelerometri 
nella prova alle vibrazioni ambientali e sono stati esportati gli spostamenti in entrambe le direzioni, 


da confrontare poi con quelli ricavati sperimentalmente, tramite la MAC. 


1° 3.05 6% 0% 3% 0% 5% 
2° 3.33 0% 0% 1% 14% 1% 
3° 4.08 3% 67% 1% 6% 
8° 5.84 1% 1% % 1% 
9° 6.28 12% 13% 1% 1% 22% 
FIHz] 2.13 2.62 3.33 3.96 4.36 
1° 2° 3° I° 6° 


Tabella 56: MAC MNI 
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4.6 Model Updating 


La procedura di aggiornamento manuale utilizzata per eseguire una calibrazione del modello è 


schematizzata in figura. 


Modello Modello 
Numerico Sperimentale 


Caratteristiche 
dinamiche numeriche 


Caratteristiche 
dinamiche sperimentali 


Modello FEM 
calibrato 


Differenza 
max. 5% 


NO 
Updating manuale 
del modello 


Selezione dei 
parametri incerti 


Modifica dei parametri 


per tentativi 


NO 


Differenza 
max. 5% 


SI 


Modello FEM 
calibrato 


Figura 152: Flowchart della procedura di aggiornamento manuale 


L’updating del modello è diviso in tre fasi: 


- Nella fase 1 viene studiato il modello MN, andando ad abbassare di volta in volta il valore del 
modulo elastico dei materiali; 

- Nella fase2 viene considerato l’effetto della rigidezza della copertura (Modello MNC) e viene 
diminuito il modulo elastico allo stesso modo della fase 1; 

- Nella fase 3 viene aggiornato il modello MNC3, inserendo le travi in legno al posto degli 


elementi shell dei solai, per andare a valutare se la soletta è collaborante o meno. 
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Alla fine, verranno messi a confronto i risultati delle tre fasi per valutare la migliore strada da seguire 
in futuro per completare la procedura di aggiornamento del modello, in modo da poter essere 


affidabile e poter essere utilizzato a supporto del monitoraggio previsto. 


4.6.1 Fase l 


Siccome il modello numerico risulta più rigido, si è scelto di svolgere altre 4 analisi con altri 4 modelli 


in cui viene diminuito del 10% il valore di tutti i moduli elastici. 


Modello 1 Modello 2 Modello 3 Modello 4 Modello 5 


MISS E E*0.9 E*0.8 E*0.7 E*0.6 
[N/mm?] [N/mm?] [N/mm?] [N/mm?] [N/mm?] 

Muratura 1 1092 982.8 873.6 786.2 707.62 
Muratura 2 1230 1107 996.3 896.7 807 

Muratura 3 1200 1080 972 874.8 787.32 

Muratura 4 1365 1228.5 LL0S.7 995.09 895.58 

Parti ripristinate 1500 1350 1215 1093.5 984.15 

Solai in legno 600 540 486 437.4 393.66 


Tabella 57: Moduli elastici utilizzati per l'updating 


Dai risultati riportati, si può notare che la riduzione del modulo elastico ha abbassato le frequenze, 
rendendo il modello meno rigido, arrivando nel modello MNS alla stessa frequenza dello sperimentale 
nel 3° modo e ad una differenza del 2% nel 2° modo. 

Dalla valutazione dell’indice MAC, si può notare un miglioramento della corrispondenza nel modello 
2, poiché negli altri modelli il valore dell’indice è inferiore. Questo ci fa intendere che forse 1 valori 
dei moduli elastici potrebbero essere corretti nel modello MN2, anche se bisognerà fare una 
valutazione precisa per ognuno di loro, tramite diverse prove, per trovare il range di valori ottimali. 
Ovvero, non basterà più diminuirli tutti dello stesso valore, ma fare delle prove variando ogni volta 
ciascun valore in modo diverso e si potrà fare sviluppando una procedura di aggiornamento 
automatico. Quindi, dai risultati ottenuti si può dire che la migliore corrispondenza tra i dati 
provenienti dal modello sperimentale e numerico si trova nel modello MN2. Infatti, in questo l’indice 
MAC risulta il migliore, con la corrispondenza del 76% del 1° modo, del 73% del 2° modo e del 72% 
del 3° modo. 


Di seguito vengono riportati i risultati di tutti i modelli in termini di frequenze e di MAC. 
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Modello 2 (MN2) 


Modo fws [Hz] fwn2[Hz] Afi [%] 
1 2.13 2.91 -27 
Z 2.62 3,21 -18 
3 3.33 3.5 -5 


Tabella 58: Confronto delle frequenze tra il modello sperimentale (MS) e il modello numerico (MN2) 


291 6% 0% 3% 
2° 321 0% % 2% 
®. 35 17% 5% 
f [Hz] 2.13 2.62 3.33 
sh pe 3° 
Tabella 59: MAC MN2 
Modello 3 (MN3) 
Modo fs [Hz] fwn3[Hz] Afi [%] 
1 2.13 2.74 -22 
2 2.62 3 -12 
3 d.35 3.69 -10 


Tabella 60: Confronto delle frequenze tra il modello sperimentale (MS) e il modello numerico (MN3) 


1° 2.74 ILVA 

3 72% 

3° 3.69 24% 3% 

fIHZ] 2.13 2.62 3.33 
1° 2° 3° 
Tabella 61: MAC MN3 
Modello 4 (MN4) 
Modo fs [HZ] fwna[Hz] Mi [%] 

1 2.13 2.6 -18 
Di 2.62 2.84 -8 
3 3.33 3.5 5 


Tabella 62: Confronto delle frequenze tra il modello sperimentale (MS) e il modello numerico (MN4) 
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1° 2.6 
2° 2.84 
3° 3.5 25% 2% 
f [Hz] 2.13 2.62 3.33 
1° 2° 3° 
Tabella 63: MAC MN4 
Modello 5 (MN5) 
Modo fws [Hz] fwns[Hz] Afi [%] 
1 2.13 2.46 -14 
2 2.62 2.69 So; 
3 3.33 3.32 0 


Tabella 64: Confronto delle frequenze tra il modello sperimentale (MS) e il modello numerico (MN5) 


1° 
2° 
3° 


2.46 

2.69 

3.32 25% 2% 

FIHZ] 2.13 2.62 3.33 


1° 2° 3° 
Tabella 65: MAC MN4 


112 


4.6.2 Fase 2 


A partire dal modello MNI, è stato valutato l’effetto della rigidezza in copertura inserendo degli 
elementi rigidi, con il rilascio alle estremità. Il modello risultante viene chiamato Modello Numerico 
C1 (MNCI). 

È stata eseguita l’analisi modale, considerando sempre i carichi da copertura, e si può notare che il 
modello risulta più rigido dei dati sperimentali e la differenza delle frequenze è maggiore nel 2° modo 
di vibrare. Procedendo con la diminuzione dei moduli elastici, allo stesso modo della fase 1, si può 
notare una diminuzione delle frequenze, mantenendo sempre la differenza maggiore nel 2° modo, e 
notando una sua distanza sempre maggiore con la diminuzione del modulo. Nel modello MNC5, 
ovvero nel modello con il modulo elastico più basso, una buona corrispondenza si ha solamente nel 
3° modo, dove la differenza è dell’1%. In termini di MAC, ci sono due modelli con i valori migliori, 
ovvero MNC2 e MNC3. Nel primo vi è una corrispondenza del 79%, 95% e 65% dei tre modi, mentre 
nel secondo i valori sono 76%, 96% e 68%. 

L'inserimento degli elementi rigidi ha comportato un miglioramento dell’indice MAC, poiché nel 
primo modello le pareti del tiburio si muovevano in direzioni opposte, mentre nel modello MNC sono 


vincolate nel movimento dagli elementi rigidi. 


Modello C1 (MNCI1) 


1° Modo 2° Modo 
Trasl. X Trasl. Y 
F3.04 Hz F3.81 Hz 


Figura 153: Forme modali e frequenze dei cinque modi di vibrare del modello sperimentale MNC 
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Modo fis [Hz] fwnic[Hz] Afi [%] 


1 2.13 3.04 -30 
2 2.62 3.81 -31 
3 3.33 4.06 -18 


Tabella 66: Confronto delle frequenze tra il modello sperimentale (MS) e il modello numerico (MNCI) 


1° 3.04 
2° 3.81 
3° 4.06 26% 0% 
f [Hz] 2.13 2.62 3.33 
1° 2° 3° 
Tabella 67: MAC MNCI 
Modello 2C (MNC2) 
Modo fs [Hz] Swn2c[Hz] Afi [%] 
1 2.13 2.86 -26 
2 2.62 3.63 -28 
3 3.33 3.82 -13 


Tabella 68: Confronto delle frequenze tra il modello sperimentale (MS) e il modello numerico (MNC2) 


1° 2.86 

2° 3.63 

3° 3.82 26% 2% 

f [Hz] 2.13 2.62 3.33 
1° 2° 3° 
Tabella 69: MAC MNC2 
Modello 3C (MNC3) 
Modo fws [Hz] fwnsc[Hz] Afi [%] 

1 2.13 2.73 -22 
2: 2.62 3.49 -25 
3 3.33 3.67 -9 


Tabella 70: Confronto delle frequenze tra il modello sperimentale (MS) e il modello numerico (MNC3) 


1 59 76% 
2° 349 96% 
3° 3.67 28% 0% 
fIHZ] 2.13 2.62 3.33 


1° 2° 3° 
Tabella 71: MAC MNC3 
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Modello 4C (MNC4) 
Afi [%] 


Modo fis [HZ] fwnac[Hz] 
Il 2.13 2.59 -18 
Z 2.62 3,399 -21 
6) 3,93 3.48 -4 


Tabella 72: Confronto delle frequenze tra il modello sperimentale (MS) e il modello numerico (MNCA4) 


1° 2.59 
2° 3.33 
3° 3.48 25% 2% 
f[Hz] 2.13 2.62 3.33 
1° 2° 3° 
Tabella 73: MAC MNC4 
Modello 5C (MNC5) 
Modo fs [Hz] fwnsc[Hz] Afi [%] 
1 2.13 2.46 -14 
2 2.62 3.17 -17 
3 3.33 3.31 1 


Tabella 74: Confronto delle frequenze tra il modello sperimentale (MS) e il modello numerico (MNC5) 


1° 2.46 
2° 3,17 
> dal 26% 0% 
FIHz] 2.13 2.62 3.33 
1° 2° 3° 


Tabella 75: MAC MNC5 
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4.6.3 Fase 3 


Siccome il modello MNC3 sembra essere il migliore, nella fase 3 si è partiti da questo modello, 


inserendo le travi in legno del solaio, per andare a valutare se la soletta è collaborante o meno. 


1° Modo 2° Modo 2° Modo 
Trasl. X Trasl. Y Tors. Z 
Fr 2.68 Hz F 3.21 Hz F 3.51 Hz 


Figura 154: Forme modali e frequenze dei cinque modi di vibrare del modello sperimentale MNT 


Modo fis [Hz] fwnt[Hz] Afi [%] 
1 2.13 2.68 -21 
2 2.62 3.21 -18 
3 3.33 Sal -5 


Tabella 76: MAC MNT 


1° 2.68 
2° 321 
3° 3.51 
SIHz] 2.13 2.62 3.33 


TE g° 3° 


Tabella 76: Confronto delle frequenze tra il modello sperimentale (MS) e il modello numerico (MNT) 


Dai risultati dell’analisi modale (Fig.154), si può notare come i valori della MAC siano peggiorati, 
ad eccezione del 2° modo di vibrare che presenta ancora una corrispondenza del 95%. Il 3° modo è 
molto lontano dai dati sperimentali, quindi la soletta è collaborante e incide maggiormente sul modo 


di vibrare torsionale. Per questo si è scelto di non continuare a svolgere analisi su questo modello. 
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4.6.4 Confronto dei risultati 


Dal confronto delle analisi nelle tre fasi, si può notare come 1 risultati migliori si abbiano nel modello 
MNC3 ovvero il modello con la soletta collaborante e le travi di copertura, con un valore dell’indice 
MAC che si aggira attorno a 95-96%. Negli altri due modi i valori dei due modelli sono simili ed 
entrambi presentano un accoppiamento del 1° e del 3° modo intorno al 24-26%. Il valore del modulo 
elastico ottimale, relativo al modello 3, è inferiore del 20% rispetto ai valori in normativa. Questo è 
dovuto alle fessurazioni diffuse presenti nella struttura. In termini di frequenza, i risultati del modello 
MNC non si discostano moltissimo dai risultati del modello MN, ad eccezione del 2° modo. Infatti, 
nel Modello MNC la differenza delle frequenze del 2° modo risulta essere sempre maggiore. Negli 
altri due modi, invece, le frequenze non si discostano molto, anzi il più delle volte assumono lo stesso 
valore. Inoltre, per entrambi i modelli si nota una progressiva diminuzione delle frequenze 
all’aumentare del modulo elastico, come si può notare in figura 155. 

Il modello MNT invece, presenta valori bassi dell’indice MAC, ad eccezione del 2° modo e quindi si 
può affermare che la soletta sia collaborante. 

Quindi, anche se il modello MNC3 risulta essere il migliore, in termini di frequenze sono ancora un 
po' lontani e quindi bisognerebbe incentrare l’analisi in un range di valori ottimali del modulo elastico, 


tramite una procedura di aggiornamento automatico. 


Confronto Frequenze MN e MS Confronto Frequenze MNC e MS 

4,50 4,50 

4,00 4,00 

3,50 3,50 
"N 3,00 N 3,00 
& CS 
2,50 q 2,50 
[=d e 
S 2,00 5 2,00 
o o 
L 150 L 
w d 1,50 

1,00 1,00 

0,50 0,50 

0,00 0,00 

0 1 2 3 s 0 1 2 
Modi Modi 
== MS —@MNI MN2 —@®—MN3 —@—MN4 —@—MNS —®-EM —@-NM1+C NM2+C —@—NM3+C —@—NM4+C  —@—NMS5t+C 


Figura 155: Confronto delle frequenze in fase 1 (MN) e fase 2 (MNC) 
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5. Capitolo 5 


Conclusioni e prospettive 


In questo lavoro di tesi è stato creato un modello agli elementi finiti della chiesa di Santa Maria in 
Via a Camerino, con lo scopo di poter essere a supporto di un monitoraggio previsto nei prossimi 
mesi. Prima di poter svolgere la sua funzione è necessario che rispecchi fedelmente il comportamento 
dinamico della struttura alle vibrazioni ambientali e quindi deve essere ben calibrato, attraverso delle 
analisi modali in cui vengono variati i parametri che la caratterizzano. Dopo aver eseguito le analisi 
di tutti i modelli (MN, MNC e MNT), si può affermare che il comportamento dinamico è influenzato 
dalla rigidezza assiale della copertura, poiché confrontando i dati numerici, il modello MNC presenta 
una MAC con indici buoni, soprattutto quello del 2° modo (95-96%), e che la soletta del solaio è 
collaborante, visto le pessime corrispondenze nel modello MNT. Inoltre, si può notare come il valore 
del modulo elastico ottimale, quello utilizzato nel modello MNC3, sia inferiore a quello presente in 
normativa, utilizzato nel modello MNI1. Questo è dovuto al fatto che nella struttura è presente uno 
stato fessurativo diffuso. 

Anche se 1 risultati della MAC sono buoni, non sono definitivi visto che un modello per essere 
calibrato al meglio dovrebbe avere una differenza massima del 5% con i dati sperimentali, per tutti i 
modi considerati. 

Il lavoro futuro consisterà nello sviluppare una procedura di aggiornamento automatico, in modo da 
rendere l’updating più veloce e preciso, con l’utilizzo ad esempio di algoritmi machine learning e 
svolgendo analisi tipo Montecarlo. In questo modo si potrà assegnare un range di valori dei parametri 
da assegnare, nel quale verranno svolte le analisi per trovare i valori corretti dei diversi materiali 


costituenti la struttura. 
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